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Wind-Tunnel Testing Problems in 
Superaerodynamics 


HSUE-SHEN TSIEN* 
Massachusetts Institute of Technology 


SUMMARY 


The problems in the experimentation of rarefied gas are dis- 
cussed. First, the extremely large viscous effects in a wind-tun- 
nel nozzle are shown. Then the difficulties of flow measurement 
are surveyed, pointing out particularly the unconventional be- 
havior of the Pitot tube in rarefied gas. The performance of a 
hot-wire anemometer is then studied in some detail to show its 
feasibility. Finally, the rules for achieving complete flow simi- 
larity of rarefied gas flow are formulated. 


INTRODUCTION 


IND TUNNELS ARE, perhaps, the most useful tool 
W:: aerodynamic investigations and certainly have 
contributed much in the modern development of fluid 
mechanics. Itis thus natural, when one turns to a new 
field of aerodynamics, the aerodynamics of rarefied 
gases or superaerodynamics, that one should think of 
using the wind tunnel again. Only here it has to be 
adapted to entirely new circumstances, and many new 
problems, both in its design and in its operation, 
appear. It is the purpose of this paper to discuss some 
of these problems so as to gain an orientation in the new 
field of experimentation. 


(1) TUNNEL DESIGN 


To test models in the wind tunnel at its test section, 
it is of primary importance to obtain a uniform stream 
at the desired temperature, pressure, and velocity. 
For subsonic wind tunnels with ordinary pressures, this 
can be achieved without much difficulty. For super- 
sonic wind tunnels at ordinary pressure, the expansion 
part of the tunnel before the test section or the nozzle is 


Presented at the Session on Aerodynamics of the Upper 
Atmosphere, Sixteenth Annual Meeting, I.A.S., New York, Janu- 
ary 26-29, 1948. 

* Professor of Aerodynamics, Guggenheim Aeronautical Lab- 
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first designed to obtain a uniform stream at its exit 
without considering the effects of the viscosity of air. 
Then the boundary layer along the wall of the nozzle is 
calculated with the pressure gradient thus determined. 
Finally, the displacement thickness of the boundary 
layer, or the needed space for the boundary layer flow 
at lower velocity, is provided by making the nozzle 
larger than the dimensions first determined by the 
calculated amount.’ This design procedure is found 
to give satisfactory nozzles for supersonic wind tunnels. 

However when one tries to apply the same design 
procedure to the superaerodynamics wind tunnel, one is 
immediately confronted with the difficulty of extremely 
large viscous effect. In other words, the boundary 
layer will be so thick as to occupy the main portion of 
the nozzle passage. To demonstrate this effect, let the 
length of the test section be Z and the width of the 
square test section be b. Then the Reynolds Number 
based upon the conditions in the test section is Re = 
UL/v, where U is the velocity in the test section. If, 
as a rough estimate, we take the thickness of the boun- 
dary layer to be zero at the beginning of the test section 
and equal to a value 6 calculated by the well-known 


. Blasius formula for a flat plate at the end of the test 


section, then 
5 = 3.65L/ VRe 


Now, if this boundary layer actually occupies half the 
tunnel width 5/2, then 


= = (2) 


On the other hand, the ratio of the mean free path / and 
the boundary-layer thickness 6 is known? to be 


1/6 = (1.255 WV 7/3.65)(M/ V Re) (3) 


where y is the ratio of specific heats and can be taken as 
1.4 and M is the Mach Number in the test section. By 


3. 3.65L/ V Re 
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Fic. 1. /, mean free path; 6, boundary layer thickness at end 


of the test section of width 6 and length L; and M, Mach Number 
of the test section. 


combining Eqs. (2) and (3), we have 
(l/5)(L/b) = 0.0557M (4) 


This relation is shown in Fig. 1. Thus for a Mach 
Number M equal to 2 and L/b = 2, the boundary layer 
will completely fill up the test section, if the mean free 
path is equal to 5.6 per cent of the boundary-layer thick- 
ness or 2.8 per cent of the tunnel width. This means 
that the extremely strong viscous effect at low densities 
makes the ordinary concept of designing a wind tunnel 
totally inapplicable. 

The extremely thick boundary layer where the veloc- 
ity increases from a small value near the wall to some 
supersonic velocity at the center of the nozzle, also 
gives subsonic velocities in a rather large portion of the 
nozzle. Since pressure disturbances downstream can be 
transmitted upstream in subsonic flows, the flow in the 
test section of a low-pressure tunnel will be sensitive to 
changes in the diffuser even if the main stream velocity 
at the center of the nozzle is supersonic. This is, of 
course, a new phenomenon in the superaerodynamics 
tunnel not found in conventional supersonic wind tun- 
nels. 

The large viscous effects can also be demonstrated by 
calculating the ratio of frictional loss on the walls of the 
test section and the shock loss in the diffuser after the ° 
test section. Consider the diffuser to be a straight tube 
of approximately the same cross-sectional area as the 
test section, then the pressure loss due to friction, Ap, is 


frictional force/b? 
(pU?/2)4bLC,(1/b?) 


Taking C; to be the Blasius value or C, = 1.328/ Vv Re, 
we have 


Api 


Ap: = 2pU%(L/b)(1.328/ V Re) (5) 


Now, if p is the static pressure in the test section, then 
the pressure for ideal isentropic compression in the 
diffuser is 
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p{1 + — - > 


If the actual pressure rise in the diffuser is estimateq as 
that due to a normal shock without further recovery 
then the actual pressure rise is 


{[2y/(y + 1)]M? — + 


Therefore, the pressure loss due to shock Ap, is 


/ 
27 


6) 


By combining Eqs. (5) and (6), the ratio of these two 
pressure losses is 


y-1 
9 


2yM2(L/b)(1.328/ V Re) 


(: 4 we)” = ( 2Qy 
(7) 


Introducing the mean free path ratio given by Eq. (3), 
we have 


y¥~-1 


Api 
Ape 
‘) 6.528yM 
(; ( (: 
2 y+1 y+1 
(8) 


This relation is plotted in Fig. 2. Therefore, if the 
Mach Number M is 2 and L/b = 2as before, then, when 
the ratio (//5) is 0.056, the ratio of frictional loss to 
shock loss is 0.628. Hence, the frictional loss and the 
shock loss are of the same order of magnitude. 
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Fic. 2. 1, mean free path; 5, boundary layer thickness at ea! 
of the test section of width 6 and length Z; and M, Mach Nur 
ber of the test section. Api, pressure loss by friction; Aps, P® 
sure loss by shock. : 
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WIND-TUNNEL TESTING 


These large viscous effects are fully confirmed by the 
recent tests on the 1 by 1-in. low-pressure wind tunnel 
at the University of California.* The test nozzle 
(Fig. 3) was designed for Mach Number 4 without con- 
sidering the viscous effect of the medium. During 
test, the static pressure on the wall at the exit of the 
nozzleismeasured. This pressure is equal to 175 and 68 
micronst for the two tests presented in Figs. 4 and 5. 
The apparent Mach Number V/ is the Mach Number 
calculated from the dynamic pressure measured by a 
Pitot tube by using the Rayleigh formula. Since there 
is the complication of large viscous effect in the Pitot 
tube reading as shown in the following section, this 
apparent Mach Number is only qualitative and cannot 
be taken as the exact value. However it is apparent 
from Figs. 4 and 5 that the boundary layer in the test 
section is indeed extremely thick and fills up the whole 
space. This large boundary-layer thickness makes the 
space available for the expansion of the central potential 
fow very small, if it exists at all. Therefore, the 
maximum Mach Number reached at the center of the 
nozzle is much smaller than the design Mach Number of 
4, At the lower pressure, the influence of slip at the 
wallisalsoevident. This has the tendency to make the 
flow more uniform. However, the extremely low 
Mach Number at the test section indicates again the 
strong viscous effect in converting much of the pressure 
energy into heat energy. 

These elementary calculations and preliminary test 
results make it clear that, for the design of the nozzle 
and test section for a superaerodynamics wind tunnel, 
it is no longer possible to separate the compressibility 


* Experimental work done under contract with the Office of 
Naval Research. The author is deeply indebted to Profs. R. G. 
Folsom and E. D. Kane for permission to use their unpublished 
results. 

t 1,000 microns = 1 mm. Hg; one atmosphere = 0.760 X 105 
microns. 


Fic. 3. Test section of University of California Low Pressure 
Wind Tunnel No. 2. 


IN SUPERAERODYNAMICS 575 


STATIC PRESSURE 
° 


68M 
T T | T | ] ] 
| [| | | 
25 ; t + 
i | | | | 
| | | | | | 
‘ 
| | |_| | 


VERTICAL POS/T/ON 
INCHES 


Fic. 4. Vertical velocity distribution in test section of Uni- 
versity of California Low Pressure Wind Tunnel No. 2. 
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Fic. 5. Horizontal velocity distribution in test section of Uni- 
versity of California Low Pressure Wind Tunnel No. 2. 


effect and the viscous effect. In fact, the concept ‘of 
boundary layer is also of doubtful value because of the 
extremely small Reynolds Number encountered. 
Therefore, to design such a nozzle to obtain the nearest 
approximation to the ideal uniform flow, it will be 
necessary to use the exact Navier-Stokes equations in- 
stead of the approximate boundary-layer equations. 
Of course, it may be argued that for superaerodynamics 
the Navier-Stokes equations are no longer exact and 
that additional corrections must be added.2 However, 
recent investigations by Schamberg* have shown that 
these additional corrections are small in the case of the 
slip flows concerned here and will not essentially alter 
the flow pattern. Hence, for a first approximation just 
like the nonviscous isentropic flow is a first approxima- 
tion for ordinary supersonic nozzles, we can use the 
Navier-Stokes equations. The simplest case to be 
considered is certainly the axially symmetric nozzles. 
If x is the coordinate in the axial direction, r is the co- 
ordinate in the radial direction, and u and v are the 
corresponding velocity components, the equations are: 


(Opu/Ox) + (1/r)(0/Or)(prv) = 0 (9) 
p(Du/Dt) = —(0p/dx) + (Grad 7), (10) 
p(Dv/Dt) = —(dp/d,) + (Grad 7), (11) 
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o(Grad 7),} + 2 ll > 7 (12) 


+ oT) = — {u(Grad 7), + 


ox or 
where 
D/Dt = u(0/dx) + v(0/dr) 
® = dissipation function 


stress tensor 


Eggs. (9), (10), (11), and (12), together with the equation 
of state, 


b/p = RT (13) 


then determine the five unknowns uw, v, p, p, and T. 
Of course, the actual process of making this calculation 
will be extremely tedious, and some approximate 
method of solution may have to be developed. One 
possibility would be to adapt the K4rm4n-Polhausen 
method for boundary layer to this case: We integrate 
the differential equations once with respect to r and thus 
only try to satisfy the equations ‘‘on the average” over 
the cross section of the nozzle. The “distribution” of u, 
v over the cross section will then be set in the form of a 
polynomial in r. Initial study by this procedure has 
already been made by Schaaf‘ at the suggestion of the 
author. 

For ordinary supersonic diffusers, high efficiency of 
pressure recovery can generally be achieved by using a 
long diffuser. However, for a superaerodynamic wind 
tannel, because of the extremely large loss through fric- 
tion, long diffusers are undesirable. In fact, the pres- 
sure loss can be reduced by using the shortest possible 
diffuser. 


(2) FLow MEASUREMENT 


The quantities that determine the flow field are two 
of the three variables p, p, and T and the velocity com- 
ponents. The quantities p, p, T are related by the 
_ equation of state, and therefore only two are necessary 
for the determination of all three. Generally, for 
wind-tunnel work, the quantities actually measured are 
p, p, and g, the magnitude of the velocity. 

For the measurement of pressure, a manometer is 
used. For ordinary pressures, one uses a fluid man- 
ometer filled with water, alcohol, or mercury. However 
for the extremely low pressure encountered in the super- 
aerodynamic flow, some other form of manometer is 
necessary. One of the most successful types is the 
Pirani gage. The conventional form of Pirani gage has 
a pressure sensitivity of about 10-? micron. It utilizes 
the change of temperature of a wire heated with con- 
stant energy caused by a change in the pressure of the 


measured by the change in the resistance of the wire. 
The wire is located in a small chamber which is con- 
nected to the point of measurement by a hole, flush with 
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gas surrounding it. The temperature change is . 
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the gas stream if static pressure is to be measured. The 
question of best design of the connecting tube for quick 
response has been studied by Schaaf. 

To measure the density p, the conventional method 
utilizes the difference in the velocity of light rays ip 
mediums of different density. With different optical 
arrangements, we have the shadowgraph method, the 
schlieren method, and the interferometer methog. 
However, if the density of the medium is low as in the 
case of superaerodynamic flow, the sensitivity of these 
methods becomes extremely poor. For instance in the 
case of the schlieren method, the percentage change in 
illumination J on passing through a region of thickness } 
is given by 
p An 
where po is the air density at 32°F. and 1 atmosphere 
pressure and where Ap/ An is the density gradient nor. 
mal to the light ray. f and e are the focal length and 
the normal unobscured width of the light source image 
perpendicular to the knife edge. & is a factor of order 1, 
determined by the particular optical path used. There. 
fore the sensitivity of the schlieren method decreases in 
proportion to the factor (p/p). Some improvement 
can be made by altering the quantities f and e, but 
practical limitations and diffraction difficulties do not 
allow the increase of sensitivity to satisfactory values, 

A new approach to this problem of density measure- 
ment is the method of absorption. It is found for in- 
stance, that oxygen at low pressures shows a strong 
absorption band at wave lengths around 1,470A. or 
ultraviolet light. The percentage absorption is pro- 
portional to the number of molecules that meet the 
light ray and is, therefore, proportional to the density of 
the gas. The measurement is then similar to that of 
the interferometer method where the density is deter- 
mined. A similar method is the utilization of the after- 
glow of nitrogen.. These methods are now being 
studied by Evans.® 

The conventional method for the measurement of 
velocity is through the use of dynamic pressure rise ina 
Pitot tube. A straightforward application of this 
method is, however, difficult for rarefied gases, since the 
formula used is based upon the neglect of viscosity 
effects. But for rarefied gases the viscosity effect is of 
great importance, as pointed out in the previous section. 
Then the dynamic pressure would be different from that 
given by the usual formula. To estimate this effect, 
let us consider the case of low Mach Number so that 
compressibility effects can be neglected. Then, asa 
first approximation, take the flow field around the Pitot 
tube as that of a source of strength S in nonviscous flow 
of uniform velocity U. The “radius” of the tube ais 


= VS/rU 


and the stagnation point is located at 
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R = VS/4xU = a/2 (15) 
The velocity introduced by the source is then 
— U(1/4)(a?/R?) 


By calculating the viscous stress from this approximate 
disturbance velocity, we have, for flow along the axis, 


+-— = pla (16) 


Hence, if po is the stagnation pressure and p° is the 
static pressure, 


R 


uf 


1 8 v 
~ pU? {1 — ——- — 17 


For rarefied gases, the value of v/aU or the reciprocal of 
the Reynolds Number of the Pitot tube could be of the 
order of unity. Then the dynamic pressure rise po — p° 
is not the usual value '/.9U? but a value much less than 
that. In fact, previous investigations by Barker’ and 
Homann’ indicate that the Reynolds Number a’U/y, 
where a’ is the radius of the mouth of tube, must exceed 
30 in order to reach the usual dynamic pressure rise 
U*. 


p= 


TABLE 1 

The Functions F, and F, [cf. Eg. (21)] 
U/ce F, F, 
0 1.77245 0 
0.2 1.73751 0.07020 
0.4 1.63880 0.27269 
0.6 1.49248 0. 58560 
0.8 1.32021 0.97843 
1.0 1.14328 1.42053 
1.2 0.97825 1.88555 
1.4 0. 83480 " 2.35492 
1.6 0.71628 2.81812 
1.8 0.62153 3.27117 
2.0 0. 54683 3.71356 
2.2 0.48790 4.14672 
2.4 0.44090 4.57300 
2.6 0.40268 4.99395 
2.8 0.37100 5.41117 
3.0 0.34420 5.87498 


When the velocity of flow is high, we have an added 
complication due to the shock. The conventional Ray- 
leigh formula for Pitot tubes in supersonic flow is based 
upon the assumption of a thin shock wave ahead of the 
Pitot tube. Now the thickness of the shock is propor- 
tional to the mean free path. Hence, in rarefied flows, 
the thickness of the shock will be so increased as to 
cause interference with flow in the neighborhood of the 
Pitot tube. This, together with the viscous effect men- 
tioned in the previous paragraph, definitely shows the 
inapplicability of the Rayleigh formula for supersonic 
velocity of rarefied gases. 


(3) Hot-Wire ANEMOMETER 


With the great complications in applying the conventional velocity-measuring device to superaerodynamic flows, 
one is naturally led to the thought of other avenues of approach. One possibility is the use of hot wires. If the 


wire diameter is of the order of 0.0001 in. and if the pressure of the gas stream is approximately 100 microns, the 
ratio of the mean free path to the wire diameter will be approximately 180. Therefore, the flow around the wire is 
definitely the free molecule flow.?. Thus we have a simple physical situation, which is an improvement over the 
rather uncertain circumstances of mixed dynamic and viscous effects for the measurement of velocity by a Pitot 
tube. It therefore seemis worth while to explore this possibility by a trial calculation of the performance of such a 
hot-wire anemometer. 

If @ is the inclination of the solid surface to a gas stream that has a macroscopic velocity U and a Maxwellian 
molecular velocity distribution, the translational energy of molecules E;, incident upon the unit area is 


Ew = Wi (« + + sin 6 c+ + erf sin (18) 


- T 


where c? = 2RT, T is the temperature of the gas stream, and erf is the error function. Now let r be the radius of 
the hot wire. Then the total energy £; incident upon a unit length of the wire is the sum of translational energy 
and internal energy. If c, is the specific heat at constant volume, this total energy per unit length of wire is 


vr \L2 2 


vi[5u + o[ 1 + sino) | ao (19) 
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The integrals in Eq. (19) can be expressed in terms of tabulated functions (see Appendix). Thus, 


E, = per HE r |r, (2) E (R+ | (20) 


c Vr —r/2 2¢ ( 1) 


F, (5) - sin + erf sin dd = A € =) (29) 


The J) and /; are the modified Bessel functions of the first kind of orders zero and one, respectively. The functions 
F, and F, are tabulated in Table 1. . 


If 7 is the wall temperature and a is the accommodation coefficient, the difference between the energy E, inci- 
dent upon the surface and the energy E, carried by the molecules re-emitted from the surface is given by 


where 


where £,, is the energy that would be carried away by the molecules if the re-emission were at the temperature 7, 
of the wire. Therefore, 


This difference of energy is then the net energy input to the wire per unit length of the wire by the air stream. 


If 7 is the electric current heating the wire and Q is the resistance of the wire per unit length at the wire tem- 
perature, the heat input per unit length of wire by the heating current is 7°2. Heat is lost from the wire by radi- 
ation. If ois Stefan-Boltzmann constant and ¢ is the emissivity of the wire surface, the radiation heat loss per unit 
length is 2rrecT.,'. Therefore, if the wire has reached a steady condition, the heat balance requires 


ap \/2RT GE U2 + Ts) | + 


Fy U?+(R+¢)(T— r.) + = QmreoT (24) 


This equation can be put-into somewhat simpler form by using the following relation: 
R= — = — 1) (25) 


Furthermore, if we take 7) as the reference temperature at which the resistance 2 is % and the corresponding 
temperature coefficient of the resistance is 8, the resistance 2 can be expressed as 


= + — To)] (26) 
Now let 
d= (7) 
Then, from Eq. (26), 
= [(\ — 1)/8To] + 1 (28) 


Now introduce pp as the reference density and % as the reference heating current; then Eq. (24) can be written as 


1+ 


The particular values of the reference temperature 7, the reference density po, and the reference current % ate 
not yet fixed. We fix these quantities now by requiring that 
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BT) = 1 
~/2e0To* = 1 (31) 
(32) 


Then Eq. (29) simplifies into 


T/T» c c y¥-1 T/T» 10 

This is then the performance equation of the hot wire in free molecular flow. 

Now let us investigate in greater detail the case of a bright platinum wire. To satisfy Eq. (30), Ty) = 492°R. 
The value for « and @ can be taken to be 0.08 and 0.90, respectively. Then Eq. (31) gives the corresponding pres- 
sure Po for po and 7% as 

po = = 3.37 microns 
Let the radius 7 of the wire be 0.0001 in. Then Eq. (32) gives the reference heating current ip as 
iy = = 0.274 milliamp. 


where the resistivity of the platinum is taken as 10.96 XK 10-* ohm-cm. Therefore the order of magnitude of the 
different quantities is entirely satisfactory. 

If the wire is used with a constant heating current, then Eq. (33) can be used to calculate the relation between the 
resistance ratio \ and the velocity ratio (U/c) at constant air-stream density and temperature.’ This is done for 
p/p = 1, T/T) = 1, and i/% = 1,* and the result is given in Fig. 6. It is seen that the sensitivity of the instru- 
ment is good. Of course, the behavior of the hot-wire anemometer will be actually determined by calibration for 
any experiment. Since the performance of the wire is strongly influenced by the accommodation coefficient a, as 
shown by Eq. (29), it will be necessary to find materials that can hold this coefficient constant for a considerable 
period of time so that no frequent calibration is required. However, the present analysis seems to indicate the 
feasibility of such an instrument for measurements in rarefied gases, and further research is definitely desirable. 


* The author is indebted to L. Mack for the numerical computations. 


(4) PARAMETERS OF FLOW functions of the molecular velocity distribution. 
Therefore the interaction of the molecules with the wall 


The two parameters that are directly connected with — . . 
th e only if the wall te ture, th tem- 
the flow field are the Reynolds Number Re, defined as 


Re = UL/v° (where v° is the kinematic viscosity, L is . 24 

the typical linear dimension of the body), and the s y 

Mach Number M° of the free stream. This is true 22 

even for slip flows and free molecule flows because of the | 

fact that the ratio of mean free path to the typical di- 20 boslincihienth VA 

mension can be also expressed in terms of the Reynolds A. A | | 

Number and the Mach Number. 18 x. | 1 
However, as the pressure or density is reduced, the VA | 

solid boundary of the flow enters actively into the flow 16 | |_| 

conditions by requiring not only that the microscopic Pa 

stream velocity be tangential to the surface but that the if VA 

interaction of the molecules and the wall be considered VA 

and that the radiation of energy to and from the wall be 12 

taken into account. The interaction of the molecules LY 

with the wall is so far expressed through the fraction s Ot 8 

of molecules that are diffusely re-emitted from the wall ufe 

and the accommodation coefficient a» It is known that Fic. 6. Resistance ratio \ of a round platinum wire of 0.0002 


both 5 : in. diameter, heated by a constant current of 0.274 milliamp. 
and @ are functions of the temperature of the in a wind stream of 3.37 static pressure and 492°R. static tem- 


wall, and there is reason to believe that they are also perature. U/c = \/y7/2 M; M = Mach Number of air stream. 
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perature, and the Mach Number of the gas above the 
wall are the same. These considerations seem to indi- 
cate then that for the model test to be similar to the 
prototype, the model must be made of same surface 
material as the prototype, the fluid must be the same, 
and, furthermore, the following parameters must be the 
same: (1) Reynolds Number, Re; (2) Mach Number, 
M°; and (3) free-stream temperature T°. 

The radiation heat loss from the surface is equal to 
eoT,' per unit area. However, if the model is sur- 
rounded by the walls of the test chamber, there is also a 
heat input due to radiation from walls of the test section 
to the model. Let us call this quantity g,.. Then the 
net heat loss per unit area of the surface of the model is 
eal! — g. This quantity can be rendered nondimen- 
sional by dividing it by p°Uc,T°. Call this new para- 
meter A,,; then 


An = [(eoTu* — gc)/p° m 


For the prototype, the heat from the walls of the test 
chamber is absent, but there may be solar radiation and 
the radiation from the earth and surrounding atmos- 
phere. Denoting this amount by g, the parameter A 
for the prototype is then 


A = (coT! — 


(34) 


(35) 


In order for the flow also to be same with respect to the 
radiation heat transfer, 


A =An (36) 


Because of the previous conditions on the Reynolds 
Number and free-stream temperature, Eq. (36) is the 
same as 


(€°T — g)/(€°T wt — Qe) = Ln/L (37) 


where L,, is the typical linear dimension of-the model 
and L is the typical linear dimension of the prototype. 
This means that the wall temperature of the test cham- 
ber must be so controlled that g, satisfies Eq. (37). 

This set of rather strict similarity rules for model 
testing in superaerodynamic flows is certainly difficult 
to satisfy. In what way the rules can be relaxed is the 
problem of future research. 


Appendix 


EVALUATION OF THE FUNCTIONS F; AND F» 


For the function 
/2 
e 


1 
= 


2 2/2) (1 20 
) (1 — cos dé 
Vx Jo 


1 —22/2 * (22/2) cos¢ ) 
¢ (A = pts 
dpe = Vr 
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where J is the modified Bessel function of first kind and 
under zero. The last step is made possible by the syb. 
stitution 26 = ¢. 

For the function F, 


/2 
F,(z) = (s sin 6) [1 + erf (z sin 6)] dé = 


/2 
I (z sin 0) erf (z sin @) dé 
—/nr2 
By the definition of the error function, we have 


/2 z sin 
F,(z) = (2/7) 2 sin 6 ds) do 


This form can be simplified by partial integration, 
Thus 


= (4/2) 22 Sor”? cos? e~* dg 
= 2? fo’? (1 + cos 26) e~* ag 
= 
2 
E (5) + cos et**/? 
2 
Therefore 


where J; is the modified Bessel function of first kind and 
order one. 
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Primary Instability of Open-Section 
Stringers Attached to Sheet 


SAMUEL LEVY* ann WILHELMINA D. KROLLt 
National Bureau of Standards 


ABSTRACT 


*A method is presented for predicting the primary instability of 
open-section stringers attached to sheet. The method extends 
the torsion bending theory to take account of the stabilizing effect 
of the sheet to which the stringer is fastened. The effects of the 
sheet in restraining twisting and warping of the stringer, as well 


as in forcing the center of twist to lie in the plane of the sheet, - 


are considered. The method is checked by comparing computed 
and observed values of maximum load for sheet-stringer panels 
reinforced by Z-section stringers. The load carried by the sheet 
in these panels was computed from a nomogram derived for a pre- 
vious investigation. Experimental data from two groups of un- 


‘published test results and from three other investigations are 


included in the comparison. 


(1) INTRODUCTION 


< ben ANALYSIS DESCRIBED in this report was con- 
ducted for the Bureau of Aeronautics, Navy 
Department, to provide a convenient method for esti- 
mating the stabilizing effect of the sheet on the buckling 
load in compression of open-section stringers in sheet- 
stringer panels and to check this:method against ob- 
served values of maximum load for such panels. 

Open-section stringers in a sheet-stringer panel 
under uniform end compressive stress may fail as fol- 
lows: 

(1) Primary instability by twisting about an axis 
parallel to the centroidal axis of the stringer without 
cross-sectional distortion. (The type of failure goes 
from combined twisting and bending to pure bending 
as the axis of twist is moved from the neighborhood of 
the stringer to infinity.) 

(2) Secondary instability by local buckling in which 
the cross section is distorted. 

(3) Primary instability by a combination of twisting 
and cross-sectional distortion. 

The mode of failure to be expected in any given case 
is that associated with the lowest computed critical 
stress. The analysis in this report is confined to rela- 
tively heavy cross sections of sufficient length to fail 
by primary instability as defined by (1). The 
possibility of modifying the analysis to take ac- 
count of primary instability as defined by (3) is dis- 
cussed. 

The theory of primary instability of columns by 
bending (i.e., twisting about an axis far from the col- 
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umn) is based on the early work of Euler and is well 
understood. The theory of primary failure by twisting 
about an axis close to the centroid of the column has 
received thorough investigation in recent years.'~* In 
references 3 and 4, particular attention has been given 
to the effect of attached sheet on the location of the 
center of twist. Additional information on the effect 
of sheet on the location of the center of twist is given 
in reference 7 for the analogous problem of bending of 
sheet-stringer panels. 


The effect of sheet on the load for primary failure of 
open-section stringers is discussed by Lundquist and 
Fligg* for the case of an I-section stringer. The re- 
sults take account of the thickness of the sheet and are - 
applicable to stringers of symmetrical cross section with 
center of twist either directly under the stringer or at 
infinity. Goodier* presents an analysis applicable to 
both symmetrical and unsymmetrical open-section 
stringers attached to sheet. His analysis takes into 
account the stiffness of the sheet in its own plane, which 
requires the center of twist to lie in the plane of the 
sheet. Weibel’ considers in: detail the restraint of 
stringers by the attached sheet in connection with the 
problem of bending of sheet-stringer panels. 


References 3 and 7 neglect the effect of buckling of 
the sheet between stringers on the primary instability 
of the stringer. This omission is probably not a source 
of serious error for the following reasons. In the case 
of thin sheet, for which sheet buckling may occur early, 
the principal effect of the sheet is considered—namely, 
to force the center of twist into the plane of the sheet. 
The theory takes account of this effect. For thick sheet, 
for which sheet buckling may either not occur at all 
or be of negligible amplitude, the bending stiffness of 
the sheet is approximately taken into account by the 
theory. 


In the present report, the theory of primary insta- 
bility of open-se¢tion stringers of arbitrary cross section 
will be extended to take account of both the bending 
stiffness of thé sheet and the axial restraint that the 
sheet provides against warping of the stringer. ° 


The condition that the center of twist must lie in the 
sheet conveniently reduces the necessary number of 
stringer section constants from the nine used by. 
Kappus? ‘to six; at the same time, three new sec- 
tion constants must be added to take account of the 
sheet. 
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Fic. 1. 


Notation for computing section constants Jr, Co, R, 


(2) TORSIONAL INSTABILITY OF OPEN-SECTION STRINGER 
ATTACHED TO SHEET 


The critical compressive stress for torsional instabil- 
ity of an open-section stringer attached to sheet will 
be computed from the equation? 


Ger = [GJr + G(Jr) sn + (Cr? E/L*) (1) 


where 
E = Young’s modulus 
G = shear modulus = E/2(1 + ») 


L = effective length of panel (one-half the actual 
length in the case of panels loaded between 
flat ends) 

v = Poisson's ratio, taken as 0.3 in this analysis 

Jr, (Jr)sn, C, Jp are section constants to be discussed. 


(2.1) Evaluation of Jr, (Jr) sr 


The constants Jr and (J7),, in Eq. (1) take account 
of the stabilizing effect due to torsional rigidity of the 
stringer and the sheet, respectively. J/7 can be calcu- 


lated by the formula 
u = uh 
Jr = S du (2) 
u= 
where 
u = coordinate along the mid-thickness of the cross 
section (shown dotted in Fig. 1); u = —u, 


at one end of the cross section, uw = wu, at 

other end of cross section, and u = 0 at rivet 
line : 

t = thickness of cross section 


If a more exact value of J7 is desired, the value obtained 

above may be corrected as outlined in Eqs. (5a) to 
_ (5d) of reference 6. Formulas for J; for Z-stringers, in 
terms of the dimensions, are given in Appendix E. 

The constant (/7),, depends on the rotation of the 
stringers relative to the sheet. In Fig. 2a, all of the 
stringers rotate in the same sense about their centers 
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of twist S, and in so doing cause the sheet to deflect 
laterally as shown by the dotted line. This type of de. 
formation is characteristic of asymmetrical stringers, 
In the case of certain symmetrical stringers of short 
length, on the other hand, the stringers tend to alter- 
nate in their direction of rotation, as indicated in Fig. 
2b. The lateral deflections, shown dotted in Fig, 2, 
vary as a sine wave along the length of the panel with 
nodal points at the ends of the effective panel length L. 

The stiffness of the sheet in its own plane forces the 
center of twist, S, Fig. 2, to lie in the original plane of 
the sheet at a distance s to the left of the rivet line. 

The number of stringers in the sheet-stringer panel 
is considered to be large enough so that the effect of 
edge conditions along the edges of the panel parallel 
to the load may be neglected. 

The contribution of the sheet to the stiffness may be 
determined approximately by calculating the amount 
of energy stored in the sheet for a given twist 8,, at the 
mid-length of the stringer using an approximate ex- 
pression for the sheet deflection. In the case of deflec- 
tions of the types shown in Fig. 2, it is shown in Ap- 
pendix A that the sheet contributes 


= Bo + (3) 
For deflections of the type shown in Fig. 2a, 


By = bh* [0.0952 + 0.289 (L/b)?] 


(3a) 


B, = 2.35 (bh*/L?) 


and, for deflections of the type shown in Fig. 2b, 


L\* b\? 
Bo = bhé | 0.288 + 0.119 (7) + 0.119 (7) ] 
bh’ LY L\* 
B, 1.175 35 175 { - 
75 + 2.35 (2) + 1.175 


(3b) 


“ks4 ~ 


(b) 


Fic. 2. Twisting of stringers and accompanying deformation 
of sheet. (a) Stringers all rotate in same direction, character- 
istic of asymmetrical stringers. (b) Stringers alternate in direc- 
tion of rotation, characteristic of some nearly symmetrical 
stringers. 
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INSTABILITY OF OPEN-SECTION STRINGERS ~~ 


where 
b = stringer spacing 
h = sheet thickness 
s = distance from center of twist S to rivet line 
of stringer (to be discussed below) 
L = effective length of panel ° 


Either Eqs. (3a) or (3b) are applicable, whichever leads 
to the lower critical stress according to Eq. (1). Eggs. 
(3a) give the lower critical stress for most asymmetrical 
(e.g., Z- and C-section) and for symmetrical stringers 
for which s and L are large compared to 6 (e.g., I- 
and U-sections). Eqs. (3b) give the lower critical 
stress for some nearly symmetrical stringers of short 


length. 


(2.2) Evaluation of C 
The term C in Eq. (1) is the ‘‘torsion-bending”’ con- 


‘stant. It corresponds to that part of the critical com- 


pressive stress caused by the resistance to nonuniform 
twist of portions of the stringer having bending stiff- 
ness. The value of C depends on the location of the 
center of twist of the stringer and upon the geometry of 
the stringer. 


A convenient method of computing C for a given axis 
of twist is presented in references 2 and 6. The method 
is used in Appendix B, together with the assumption 
that the center of twist lies in the plane of the sheet. 
Appendix B leads to the following procedure for com- 
puting C: 


(1) Draw a line along the mid-thickness of the 
stringer cross section (shown dotted in Fig. 1). Let u 
be the coordinate along this line. Let the value of u 
be zero at the rivet line and let u increase to the right 
as viewed at the rivet line. 


(2) Determine, as a function of u, the distance r 
(Fig. 1) from the rivet line a to the tangent line bc. 
The direction of the tangent line is the same as the 
direction of wu. The value of r is taken as positive when 
line ad must be rotated counterclockwise 90° to assume 
the direction of the tangent line bc; it is taken as nega- 
tive when line ad must be rotated clockwise 90° to as- 
sume the direction of the tangent line. The value of r 
shown in Fig. 1 is positive. 

(3) Compute the unit warping displacement w as a 
function of u from 


w= fi" rdu (4) 


(4) Compute the torsion-bending constant Cy for 
twisting about a, Fig. 1, from 


1 u = 2 
a= f wtdu — 


(5) 


where 
—u, = lowest value of u 
u, = highest value of u 
t = thickness of stringer 
A = cross-sectional area of stringer 
@ = cross-sectional area of sheet per stringer 


(for unbuckled sheet in a wide panel, 
a = bh; for buckled sheet in a wide panel, 
a = bh/2. b = stringer spacing; h = 
sheet thickness). 


(5) Compute the warping moment R from 
= tdu — t 
“= u, 
( f yt an) (6) 
—U, 
where y = coordinate perpendicular to sheet at rivet 


line, Fig. 1. 
(6) Compute the moment of inertia J from 


tau — yt du 
A+a — 


(7) 


‘x = coordinate parallel to sheet at rivet line. The 


x-axis is tangent to the mid-thickness line at 
u = 0. 

(7) The torsion-bending constant: C required in Eq. 
(1) is given by the equation 


C = Co + 2Rs + Is? (8) 


where s = the distance along the x-axis from the 
center of twist to the rivet line, Fig. 1. The value of s 
is positive if the center of twist is to the left of the rivet 
line as shown in Fig. 1. The determination of s will be 
discussed below. 

Formulas for Co, R, and J for Z-stringers are given, in 
terms of the dimensions, in Appendix E. 


(2.3) Evaluation of Jp 


The term J, in Eq. (1) is the polar moment of inertia 
about the center of twist. It is shown in Appendix C 
that J, may be computed from: 


Jy = Jn + 2M,s + (A + a)s? (9) 
where 
M = tx du (9a) 
Jy = polar moment of inertia of stringer cross 
section about rivet line 
t, A,a defined in connection with Eq. (5) 
) defined in connection with Eq. (8) 


Appendix E gives formulas for M,, Jyo, and A for 
Z-stringers in terms of the dimensions. 
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(2.4) Evaluation of s and of cc, 


The critical compressive stress for a given location of the center of twist is given by Eq. (1). In that equation 
(Jr)sn, C, and J, are quadratic functions of s determined by Eqs. (3), (8), and (9), from the panel dimensions; 7 
is the effective length of the panel; E is Young’s modulus; and G is taken as E/2.6. After substituting for J,, 
(Jr)sn, C, Jp and factoring out £, it is found that Eq. (1) can be expressed as 


(0.385; + 0.385By + 9.87 + (19.74 + (0.3858, + 9.87 
Jo + 2M,s + (A + a)s? 


In Appendix D, it is shown that the buckling stress o,, will be a maximum or a minimum if the distance s to the 
center of twist is a root of the equation 


cr 


(10) 


2 
(;) (A +) — M, (0.385B,L? + | + (;)[ + a) (0.385 + 0.385BoL + 9.87 


—Jyo (0.3858. + 9.87 + (0.385. + 0.385By + 9.87 — (987 =0 (11) 


The roots s of Eq. (11) can be determined algebraically. Substitution of s into Eq. (10) gives o,,/E. 


(2.5) Plastic Buckling mated from that in the elastic range (6¢,)elastic, given by 


It has been shown by Shanley,’ in the case of col- Eq. (1), by solving 

umns, that plasticity causes the effective modulus to Ger = (Cer)etastic(Er/E) (12) 

drop to a value between the tangent modulus and von- oe en , 

reduced modulus. Stowell’ shows how Se = 

the effective modulus in the case of plates depends on 
value of E,/E can be computed knowing Young's 

the stress-strain curve of the material and the geometry ‘ 

modulus and two secant yield strengths o;, o2. These 

strengths are defined as the stresses at the 

theories to take account of instability of the type con- 

intersections of the stress-strain curve and lines through 

ty the origin with slopes of 0.7F and 0.85£, respectively. 


modulus in the plastic range on some other basis. It is The computation is made with the help of the nomo- 
shown in Fig. 10, reference 6, that, for 125 unsupported a al 
st gram taken from reference 6 and shown in Fig. 3a or 
open-section columns that failed by torsional instability, : 
3b, proceeding as follows: 
buckling stresses computed with the reduced modulus : 
Draw a line from On Scale 1 through the 
agreed, within the experimental scatter, with the meas- : ; é 
‘ i value o,/o2 on scale 2. From the intersection of this 
ured buckling stresses. On the assumption that _,. 
, we i line with scale 3, draw a line through the value o;/a on 
stringers attached to sheet have a similar behavior, the 
scale 5. The intersection of the dast line with scale 4 
P is the desired value of E/E. 


sear 2 seme scars (3) COMPARISON OF COMPUTED WITH EXPERIMENTAL 
L RESULTS 


A comparison of computed critical stresses with ex- 
perimental values was made in view-of the simplifying 
assumptions made in the preceding analysis. Foremost 
among these assumptions are: (1) that the stringer is 
sturdy so that failure is by primary twisting instability 
with little or no cross-sectional distortion, and (2) that 
the effect of the bending stiffness of the sheet in tor- 
sionally restraining the stringer is adequately accounted 
for by the simplified analysis in Appendix A. 

An additional assumption, frequently made in using 
analyses of the type presented here, is (3) that the 
critical stress for primary instability and the failing 

Nemogram for determining reduced modulus stress are the same. This assumption will be essentially 
Fic. 3a. Nomogram for determining reduced moduhus, correct for Z-section stringers but may give a conserva- 
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Nomogram for determining reduced modulus. 


Fic. 3b. Nomogram for determining reduced modulus. 


tive value for other types of stringers, such as thin- 
‘walled angles, for which the critical stress can be con- 


siderably less than the maximum stress (reference 6, 
Figs. 11 and 12). 

Values of critical stress were computed for the com- 
parison not only as outlined above but also taking 
(Jr), = 0. This was done for panels with relatively 
stiff sheet to approximate the effect of a type of cross- 
sectional distortion localized in the web of the stringer 
in which the free flange of the stringer rotates relative 
to the flange attached to the sheet. For panels with 
sheet that is not so stiff, the effect of (J7),, should be 
small, and the calculations were intended to show 
whether the error due to neglecting (J7),, for such panels 
was negligible. 

The comparison was extended over a sufficient num- 
ber of references among the many available in the 
literature to cover the range of panel dimensions and 
stringer shapes most commonly used. The results of 
the comparison are given below for each reference in 
turn. 


(3.1) Unpublished Test Results—Group A 


TestS were made of 18 24S-T aluminum-alloy panels 
with Z-section stringers. The panels covered a range 
of lengths from 8 to 24 in.; stringer spacings of 2, 4, and 
8 in.; and sheet thicknesses of 0.025-, 0.051-, and 
0.081-in. The stringer cross section is shown in Fig. 
4a. The panel dimensions are given in Table 1. The 
panel constants, computed as outlined in this paper, are 
listed in Table 2. The area of sheet acting with each 
stringer, when the sheet was unbuckled, was taken as 
the total sheet area divided by the number of stringers 
and, when the sheet was buckled, as one-half this area. 
The load carried by the stringers, computed by con- 
sidering o,, to be the maximum stress of the stringers, 
was calculated from the panel constants and Eqs. (11), 
(10), and (12). Sheet loads were obtained by use of the 
nomogram of reference 10. 
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Fic. 4. Cross sections of stringers used in panels for which 


maximum loads or stresses were computed. All dimensions in 


inches. 
reference 12; (d) unpublished tests, Group B; (e) reference 13. 


(a) unpublished tests, Group A; (b) reference 11; (c) 


Table 1 - Dimensions of paste 


lished test results ~ 
[Panel] Length] Width ick- | Bo. of itringer vet | Area of] Area of] Area o 
ness etringers spacing | spacing) sheet each panel 
stringer 
(any) | in.) | (an) (432) 
2 17.08 -02 5 5 137 1.102 
3 16.0 16.83 | .0266 5 135 1.122 
20.00 16.8 5 5 139 1.114 
24.05 | 16.8 .0250 3 4 2 421 2136 /.1.099 
12.05 16.90 | 6 421 
7 124.05 16.91 | .026! 3 6 5 
|12.07 16.92 2 2 423 1. 
9 j24. 16.87 +0250 9 2 422 -138 1.667 
10 16.90 9 2 1.0 875 «137 2.105 
11 16.81 0512 5 1.0 +141 1.567 
12 ae) 16.80 +0510 3 3 1.0 7 1} 1.261 
1 16.87 +0813 9 2 1.0 1.372 2.593 
1 24. 16.87 5 1.0 1.366 137 2.050 
1 24.01 16.87 0814 1.0 3.373 1.767 
1 24.03 {12.91 | .0 a 327 862 
17 [24.03 8.7% | .02h8 3 4 5 +217 62 
124,02 4,86 | 4 25 121 213 
Reference 11 
Stringer A 
PB-15) 24 20.8 
Stringe 
7025 3 2220 7698 
| Po-12) 24 | 3 6 275 | 792 
Reference 12 
t 19 16 +070 +875 [1.12 
Unpublished test resulte -Grour 
Sheet No of | Stringer Rivet] area of] area of Aree of 
Penel| Length] Width! Thick- etringere| spacin spacing sheet each panel 
| (10,)] news (1n.) stringer 
ine) (ins) (ise) (ini) 
3 15.97 | 16,04) 0.9517 3 4.00 1,00 |0,829 1,404 
4 | 19.97 | 16.04] .os18 3 4.00 1.00 | ,631 193 1,410 
23,32 | 16,01] 0506 3 4,00 1,90 g10 195 1.295 


Reference 13 


Stringer A, Figure 4(e) 
1 23.24 26.10] 0.1016 € 5.08 1.25 | 2,652 180 3.73 
2 | 41,52 | 26,10) .1016 6 5.08 1.25 | 2,652 180 3.73 
3 | 58,26 | 26,10) .1016 6 5,08 1.75 | 2,652 180 3.73 
4_ | 60,56 | 26,10) 6 5,08 1,25 | 2,652 2180 3.73 
Stringer B, Fi;ure 4(e) 
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Table 2 - Values of constants for 2 aga 


ublished teet resuite ~ Grou 
Sheet Stiffener 

000084: | .004010 900332 02903 | .0001425 | 08587 
+00009! 200000275 | .00%011 000342 02922 2425 | 04523 | 
000117! 00000145 | .004011 000349 02934 | .0001425 | 08547 
+0001692 | .00000102 | 000325 02892 | .0001425 | 086) 
00003: joo =.00%020 000507 032 0001425 | 08687 
20001101 | .00000239 | .00%023 900553 03279 2425 | 587 
-0000852 | .0000020 004001 000175 02639 | .0001825 | 
0003294 | .00000051 -008001 000178 02643 | .0001425 | 08587 
200000852 2004018 02 0001425 04523 08587 
001852 +00000873 | .00%0; 000! 93441 | .0001825 | 
-0007938 | .00001724 | 00085 03795 | .0001485 | 08587 
01122 00001765 | .004043 | .0008 03858 | .0001425 | 085 
005755 +0000 345) 2004080 | .000843 | .03771 | .0001825 | | 
003219 ~000070 008054 +001075 204168 | .0001825 | 08587 
0001751 +004010 .000340 02919 | .0001825 | | 
+0001649 | .000 #008006 000270 425 | | 
+0001688 | .00000102 | .00%00% | .000229 02730 | .0001425 | 

Reference 11 
Reference 12 
| te osbines cone:ante | Stiffener c newrte 
T 0.002367 | 000003573 | 0.001731 [0.0012 0.02760 7.05017 
6 093076 | 001655 21565 2 
Unpudlished test results - Grow 

3 00002037 06217 

00001311 21222 

my 06217 222 


Combined constante same as for panel 4 shen (J,),, 0 a8 sheet fully effective 
“ Gomb ined constants same as for panel 8 when (Jy),, = 0 ae sheet fully effective 
Reference 13 (continued) 


Stringer C, figure 


| 0, 91053 | 9,0000045, 0,004001 | 0.003018 | 0. 05552 


3 constants sane as for panel 10 gn a8 sheet fully effective 


Comparison between the calculated and experimental 
maximum loads (Table 3) shows that the calculated 
total loads were from 5 per cent too low to 17 per cent 
too high. For those panels where the sheet thickness 
is considerably less than the stringer thickness, 0.025-in. 
sheet, 0.055-in. stringer, however, calculated and experi- 
mental loads differed by 8 per cent or less, with the 
exception of the shortest panel (8 in.) in which the cor- 
rection for reduced modulus is relatively large, leading 
to a relative large effect of errors in reduced modulus. 
For those panels for which the sheet thickness was 
approximately equal to, or greater than, the stringer 
thickness, the calculated loads were from 10 to 17 per 
cent higher than the experimental loads. These panels, 
after failure, showed that appreciable cross-sectional 
distortion accompanied primary failure. This distor- 
tion probably accounts for the lower experimentally 
observed stress at failure. 

The total maximum loads were also computed with 
the assumption (J7),, = 0. The resulting values are 
also given in Table 3. - Comparison of these values with 
those obtained when (J7);, was determined by Eq. (3) 
shows that the neglect of (J7),, changed the computed 
load less than 8 per cent for the twelve panels using 
0.025-in. sheet. For the remaining six panels with 
0.051- and 0.081-in. sheet, the neglect of (J7),, leads to 
much better agreement between calculated and experi- 
mental loads. For these six panels, appreciable cross- 
sectional distortion of the stringer accompanied failure. 
The distortion was of such a nature that the lower edge 
of the web acted like a hinge. Apparently, the neglect 


1948 


of (Jr)s, can approximately account for this type of 
distortion. 


(3.2) Reference 11 


In reference 11, results are presented for aluminum. 
alloy panels with 24S-T aluminum-alloy sheet and 
Z-, S-, C-, and U-section stringers of 24S-RT aluminum 
alloy. 

Calculated results were obtained for several panels 
with Z- and S-section stringers for comparison with ex- 
perimental results. The dimensions of these panels 
are given in Table 1; the stringer sections are shown in 
Fig. 4b. 

The panel constants for the Z- and S-stringers were 
computed by the simplified formulas in Appendix E, 
considering a to be one-half the sheet area per stringer. 
Values are given in Table 2. The computed total maxi- 
mum loads, using the nomogram of reference 10 to get 
the sheet loads, are given in Table 3. Reference 11 does 
not contain compressive stress-strain curves for the 
24S-RT stringer material which might be used for de- 
termining the reduced modulus. Compressive stress- 
strain curves for typical 24S-RT material could not be 
used, since the tensile properties of the stringer mate- 
rial were considerably lower than those for typical 
24S-RT material. The tensile properties approxi- 
mated those for 24S-T material. Hence, it was decided 
to determine the reduced modulus from compressive 
stress-strain curves of 24S-T material with the same 
tensile yield strength as the stringer material. 

Comparison of the computed and experimental 
maximum loads for these panels shows that for panel 
PB-15, which was experimentally observed to fail by 
torsional instability, the computed load was 7 per cent 
too high; for panels PB-14 and PC-12, which were ob- 
served to fail by a combination of torsional and local 
instability, the computed loads were 3 per cent high and 
6 per cent low, respectively; while, for panels PB-1 
and PC-1, which were observed to fail by local insta- 
bility, the computed loads were 10 and 2 per cent high, 
respectively. Thus it is seen that the computed loads 
for the three panels that failed entirely or partially by 
torsional instability agree with the experimental loads 
within 7 per cent. For the two panels that failed by 
local instability, the computed loads for torsional in- 
stability in one case, panel PB-1, is definitely higher 
than the observed load and, in the other case, panel 
PC-1, agrees within the possible errors. The results 
for these two panels are in agreement with the condi- 
tion that the computed load for torsional instability 
should exceed the observed load when the failure is due 
to another type of instability. 

When (Jr),, was neglected, the computed loads for 
these five panels were from 10 per cent higher to 9 per 
cent lower than the experimental values. The 
neglect of (Jr), had only a minor effect on .the 


computed loads. 


(3.3) J 

Tw 
tefere 
in thi: 
shown 
and tl 


perim« 
puting 
24S-T 
in obt 
fabric: 


hes 
- 
H 
10 
12 
ie 
1 
| 
Reference 13 
| 00003273 | .12733 | .0002455 | .05750 | .1906 
adh 02622 200002905 009621 02008 12783 0002455 | .05750 21986 
05112 00001476 | .02208 | .12733 | .0002655 | .05750 | .1986 
09 90000 77; 4282 0002455 95750 
tringer figur 
5 008282 0003737 01096 02426 13733 0002483 | ,05969 
6 02391 20001168 201096 02426 13733 0002483 | .05969 2007 
7 04711 0000578 01096 | .02426 | .12733 | | .05969 | 2007 
01355 02703 | _,15024 202483 | ,05969 2007 
’ 
: 
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Table 3 - Comparison of Computed and Experimental Results for Panels 
Unpublished test results - Group A 7 
—— | Differente 
6,100 29,200 300 5,600 33,200 38,800 9.9 5,600 33,100 38,700 9.6 
| 63200 | 28,000 38200 | 29,900 35,200 | 2.9 51300 35,300 | 3.2 
»700 | 23,700 31,400 ,000 »600 31,600 6 ,000 25,800 31,800 1.3 
H hrs 23,200 30,000 5,500 »000 29,500 - .2 5,500 23,500 29,000 | = 3-3 
5,600 22,900 28 ,500 5,500, 22,100 ° - 3.2 5,500, 20,900 26,400 | - 7. 
2 3,500 15, 19,100 3,5002 16,700 20,200 5.8 3,500 16,600 20,100 5.2 
1 3,000 14,000 17,000 3,0002 13,100 16,100 - 5.3 3,000° 12,600 15,600 | = 8.2 
8 8,700 1,400 60,100 8,400 232700 60,100 0.0 8,400 51,800 60,200 
9 »200 » 300 +500 »400 »500 50,900 2.8 8,400 38,400 46,800 | - 5-5 
10 24,800 | 39,700 64,500 | 26,700 47,200 »900 14.6 26,700 37,900 64,600 
ll 17,500 21,200 »700 17,200, 27,100 » 300 14:35 17,200 22,200 39,400 1.8 
1 10,100 12,600 22,700 10,100 15,200 25, 300 11.5 10,100 13,200 23, 300 2.6 
1 » 300 »200 91,500 »400 51, 100,800 10.2 49,400 40, 300 89,700 | = 2.0 
1 3,100 23,400 6,500 36,700, 29,300 66,000 16.8 36,700, 22,300 3300 4.6 
1 1,100 11,400 2,500 21,100 15,200 36,300 11.7 21,100 12,900 »000 4.6 
1 " 17,900 22,500 » 300 18, 300 22,600 4 » 300 17,200 21,500 | - 4.4 
17 3+J00 14,000 16,700 2, 13,800 16,600 |= .6 2,800. 13,100 15,900 | - 4.8 
3540 10,540 2,0002 11,400 £.2 2,000 2.5 
Reference 11 
Stringer A, Figure 4(b) 
PB-1 ania — 19,000 2,600 18,400 21,000 10. 2,600 18,400 21,000 20-3 
PB-14) --- 17,740 2,600 18,200 2. 2,600 15,600 18,200 2. 
PB-15|  --- aus 17,300 2,500 16,100 18,600 2,50 15,800 18,300 5.8 
Stringer B, Figure 4(b) 
a _ — 21,700 2,500 19,700 22,200 2. 2,500 19,400 21,900 0.9 
20,450 16,700 19,300 ~ 2,600 16,000 18,600 | 9.0 
ference 12 
22,400 14,100 36,500 22,4003 1 37,700 3.3 22,4003 13,200 35,600 | - 2.5 
| 14,100 17;700 | - 3.8 43300 | 12,900 17,200 | = 6.5 
14,700 | 21,200 »900 | 14,700 22,100 36,800 a. 
13,300 20,700 32300 3 33,900 3 
5 13,300 | 19,300 32,600 | 13,3003 19,000 32,300 | - .9 


1 (Total 1084) oaic, (Total load)exp. x 100 
(Total Toad), 


2 Experimental value used since panels had only 1 or 2 sheet bays which was considered ta few for the nomogram of reference 10 


3 Experimental value. 


Reference 13 


Stringer A, Figure 4(e) 


etress| Calculated stress, (JT)sh gonsigered ted stress, 
iff. in 4ff. in 
Panel] Sheet Stringer Avera,e,| Sheet, Stringer,| Avera,.e avera,e Sheet Stringer | Average } average 
pei psi psi psi stress, % pes psi pei stress, 4! 
1 aa -- 28,000 27,500 4g,400 32,900 17.5 24,500 29 ,800 25,900 = Te 
—_ 25,000 | 27,500 40,400 | 30,800 13,100 13,100 | 13,100 
FH --- -- 20,400 6,000 31,800 27,500 & 9,100 9,100 9,100 5.4 
--- --- 13,500 | 13,300 13,300 _ | 13,300 1.5 2.100 2,100 1,100 7.4 
Stringer B, Figure 4(e) 
H es — 25,900 | 28,000 48,500 | 31,800 22.6 28,000 36,300 | 29,600 14.3 
--- _— 24,000 28,000 43,700 30,900 28.7 22,500 25 ,800 23,100 - 3.8 
7 --- --- 22,200 | 28,000 36,100 291900 32.9 11,000 11,000 | 11,000 =30-9 
& --- --- 13,600 | 14,100 14,100 | 14,100 3.7 1.200 1.200 1.200 7-1 
Stringer C, Figure 4(e) 
9 | --- | --- | 34,100 35,000 | %2, 000 [38,500 | 12.9 } 19,100 19,100 “44.0 
10 --- 17,200 118,700 18,700 | 18,700 8.7 10,200 10,200 _ | 10,200 —40.7 


1 (Total load) 44, - (Total load). = 100 
(Total load) 


exp. 


2 Experimental value uses since penels had only 1 or 2 sheet baye which was considered too 


reference 10. 
3 Experimental value. 


few for the nomogram of 


(3.3) Reference 12 


Two of the 24S-T aluminum alloy panels reported in 
reference 12 were analyzed by the method presented 
in this report. The stringer section for these panels is 
shown in Fig. 4c, the dimensions are given in Table 1, 
and the panel constants are given in Table 2. The ex- 
perimentally determined sheet load was used in com- 
puting the total load for panel 1, which had Alclad 
24S-T aluminum-alloy sheet; the nomogram!® was used 
in obtaining the sheet load for panel 6, which was 
fabricated with 24S-T sheet. The experimental and 


computed values are given in Table 3. They agreed 
With (Jr), = 0, the calculated 


within 4 per cent. 


loads were 3 and 6 per cent too low for panels 1 and 
6, respectively. 


(3.4) Unpublished Test Results—Group B 


Five Alclad 24S-T aluminum-alloy panels were tested 
to study the effect of length on the ultimate load of the 
panels. The comparison of computed and experimental 
loads was carried out for the three longest ones, panels 
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3, 4, and 5. The dimensions of these panels are given 
in Table 1. - 

The stringer section for these panels is shown in 
Fig. 4d and the panel constants are given in Table 2. 
The computed maximum total loads, Table 3, were 
from 4 to 6 per cent higher than the observed loads. 
When (J7);, was neglected, the computed total loads 
differed from the observed loads by —1 to 3 per cent. 


(3.5) Reference 13 


In reference 13, test results are reported for 576 
24S-T aluminum-alloy panels with formed Z-stringers. 
The panels cover a wide range of lengths, stiffener 
spacings, and ratios of sheet thickness to stiffener 
thickness. Of these, ten panels were selected, repre- 
senting three types of stringers (see Fig. 4e). The 
dimensions, computed from the nondimensional form 
of reference 13, are given in Table 1. The thickness of 
the sheet used in these panels was equal to, or greater 
than, the stiffener thickness. The stiffeners were 
formed from sheet and had no bulb, lip, or other stiffen- 
ing element at the end of the flange. In addition, the 
outstanding flange was narrow. 

The maximum stringer and sheet stresses were com- 
puted from the preceding theory and the nomogram of 
reference 10. Tke panel constants are given in Table 
2. The median stress-strain curve for the stiffeners 
was used to correct for the reduced modulus. The 
stresses are listed in Table 3. The computed average 
stresses for the panels were 13 to 35 per cent greater 
than the measured stresses, except for the longest panels 
for which they were from 2 per cent too low to 9 per 
cent too high. The discrepancy may have been due to 
local instability of the stringers or to cross-sectional 
distortion accompanying primary failure, particularly 
for stringers A and B, which were attached to heavy 
sheet. 

With (J7),, = 0, the computed stresses were from 
55 per cent less to 14 per cent greater than the ob- 
served values. The lack of agreement in this case would 
indicate that cross-sectional distortion loca'ized in the 
web of the stringer was not predominant in these 
stringers with their narrow outstanding flanges. 


(4) ConcLusIoNns 


A comparison of computed and observed maximum 
loads for panels reinforced by Z-stringers of several de- 
signs leads to the following tentative conclusions: 

(1) For panels reinforced by Z-stringers having a 
stringer thickness of two or more times the sheet 
thickness, the proposed method of computing stringer 
loads can be expected to give total loads accurate to 
within about 8 per cent. The computed total loads 
were lowered less than 8 per cent with the simplifying 
assumption (J7), = 0. 

(2) For panels reinforced by Z-stringers having a 
bulb. or lip at the end of the outstanding flange and 
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having a stringer thickness equal to, or less than, the 
sheet thickness, the proposed method of computing 
stringer loads may give total loads up to 17 per cent 
too large; however, if used with (Jr), = 0, it can be 
expected to be accurate to within about 8 per cent in 
many cases but cannot be relied on in all cases, 
(3) For panels reinforced by Z-stringers having nar. 
" row outstanding flanges without a bulb or lip and hav- 
ing a stringer thickness equal to, or less than, the sheet 
thickness, the proposed method of computing stringer 
loads may give total loads up to 35 per cent too large, 
With (Jr), = 0, the computed loads may be as much 
as 55 per cent low and as much as 14 per cent high. 


APPENDIX A—EVALUATION OF (J7)<), 


In order to evaluate (J7),, for the case where all the 
stringers rotate in the same direction, Fig. 2a, deter- 
mine the energy stored in the sheet when the stringers 
have rotated through an angle @ about a center of twist 
S, Fig. 2a, where S is a distance s to the left of the rivet 


line. The following notation is used: 

G = shear modulus = E/2.6 

D = flexural rigidity of sheet = Eh*/10.92 

s = axial coordinate along stringer 

x = transverse coordinate in sheet, with origin at 
rivet line 

w = lateral deflection of sheet 

b = stringer spacing 

h = sheet thickness 

L = effective panel length 

6. = maximum twist of stringer relative to ends of 
panel of length L 


Approximate the lateral deflection of the sheet by 


3x? . b 
OS*S>5 
2 b 
w= + 2) 3730 


(1A) 


From page 307 of reference 14, the energy stored in the 
sheet is 


9 = \\Ox? 


: Oz? 
\2 
2A 


Substituting Eq. (1A) into Eq. (2A) and performing the 
integration gives 


= 2 — 
V = On? (s + 0.987 + 0.116 + 24.47, 
(3A) 


The energy stored in a rod of torsion constant (Jr)n 
and length L twisted through angles 


| 
= 
UJ 
al whe 
strit 
dete 
|| 
B, 
a The 
using 


1 at 


s of 


2A) 
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6 = 0» sin (r2/L) (4A) 


V! = (w?/4L)G(Sr) sn Om? (5A) 
Equating V and V’, Eqs. (8A) and (5A) and solving for 


(Jr)sn gives 


L\2 2 
Tr)sn = Oh? | 0.0952 + 0.289 + 2.35 + 
( r)s x, 


b 
0.0112 (6A 
) 


Since b will ordinarily be less than L, the last term in 
Eq. (6A) can be neglected and Eq. (6A) rewritten 


(Jr) sn = Bo + B,s? 


where 
‘By = bh*®[0.0952 + 0.289(L/b)?] 

(7A) 
. B, = 2.35bh*/L? 


In the case shown in Fig. 2b where the direction of 
stringer rotation alternates from stringer to stringer, 
determine the energy stored in the sheet when the 
lateral deflection is approximated by 


2 
w= + () sin (sin = + (8A) 
L b 


where = tan~! (s/b). 
Substituting Eq. (8A) into Eq. (2A) and performing 
the integration gives 


2 


1 1 
V=0,'D + 5?) + Dp + x) (9A) 


Equating V, Eq. (9A), to V’, Eq. (5A), and solving for 
(Jr) sn gives 


2 
(Jr)in = bh’ {0.238 + 0.119 + 0.119 (;) + 


ll 1.175 + 2.35 (7) +°1.175 (7)'}} (10A) 


Eq. (10A) can be rewritten 
(Jr)sn = Bo + B,s? (114A) 


where 


L\? b\2 
By = bhi [0.28 + 0.119 (F) + 0.119 @ | 
bhi _ (LN? __ 


APPENDIX B—EVALUATION OF C 


The torsion-bending constant C will be evaluated 
using the notation of Section (2.2) and Fig. 1. The 


(1) The sheet has sufficient shear stiffness in its 
plane to force the center of twist to lie in the plane of 
the sheet. 

(2) The axial rigidity of the sheet restrains warping 
of the stringer as though an additional area a were 
added to the stringer at the rivet line. This area a 
would equal the sheet area per stringer, bi, in the case 
of unbuckled sheet, and it would equal the effective 
incremental sheet area per stringer, bh/2, (correspond- 
ing to a small change in load) in the case of buckled 
sheet. (Fig. 6 of reference 15 gives about 0.4bh for 
simply-supported sheet, while Fig. 4 of reference 16 
gives about 0.5bh for clamped sheet. A value of 0.5dh 
was used in all cases, since this was considered a suffi- 
cient approximation.) 

Using Eq. (28) of reference 6 and the notation of the 
present report, the torsion-bending constant C about 
an arbitrary center S, Fig. 1, can be computed from the 
formula 


1 
du — 3 


(1B) 


where 
wv, = Je" r, du (2B) 


r, is the distance from the center of twist S, Fig. 1, 
to the tangent line. It should be noted that w, = 0 
at the rivet line where wu = 0. From Fig. 1, 


r, = r + s(dy/du) (3B) 
Substituting this in Eq. (2B) gives 
w,= fo'rdut So" s (dy/du) du 


giving 
WwW, = wt ys (4B) 
where 
w= du (5B) 
Substituting Eq. (4B) into Eq. (1B) gives 
C = Co + 2Rs + Is? (6B) 
where 


1 2 
= 2 


u= Uh 1 

u= UA u= Uh 
( f wt au) ( i yt au) 

—ul “= 
1 2 

= 2 


The contribution to C of the bending stiffness in a 


following assumptions are made regarding the effect of direction normal to u is unaffected by the presence of 
the sheet on the constant C: 


the sheet or the location of the center of twist. It can 
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al “tw 


Rivet \e | Rivet line 

L b,/2-] = 

“1b, /2- 


Fic. 5a. Stringers with equal flanges and lips. Fic. 5b. Stringers with unequal flanges and no lips, 


be computed from Eq. (18), reference 6. The corre- AppeNpIx D—EVALUATION OF 5 
sponding term in Eq. (1) will be small compared to the ‘ 

other terms, except possibly in the case of short panels In Eq. (10) the value of o,,/E is given as a function 
using angle stiffeners. In doubtful cases, it should be 0f the section constants and of the distance s from the 


computed and added to the value of C given in Eqs. (8) center of twist to the rivet line. Both the numerator 
and (6B). and denominator of the right side of Eq. (10) are func. 


tions of s. Let f,(s) be the numerator and f(s) be 
APPENDIX C—EVALUATION OF Jp the denominator; then 
The quantity J, is defined as the polar moment of Oer/E = fa(s)/fals) (1D) 
inertia of the section about the center of twist S, Fig. 1. 
The sheet is taken into account by adding an effective (2) fa(s) dfn(s) — fy(s) dfa (s) 
sheet area a to the stringer at the rivet line. Then E - ds ds (2D) 
ds [fa (s)]? 


Jp = *h + x)? + du + as? (1C) 
J 2M,s + (A + a)s? 9c) When o,,/E is either maximum or minimum d(¢,,/E)/- 
ds = 0. Applying this to Eq. (2D), 


Jo = — fals)[dfels)/ds] = 0 (8D) 
= ft xt du Substituting the numerator and denominator of Eq. 
“= —M (10) for f,(s) and fe(s), respectively, in Eq. (3D) gives 
du Eq. (11). 


APPENDIX E—CONSTANTS FOR PANELS WITH Z-STRINGERS 
(1) Stringers with equal flanges and lips (see Fig. 5a): 


Jr = “/sltw* (2br + bw) + (2bz)] 
A tr + tw[2bp + (bw/cos 6) 


2h 3 
[tw (bpbw? bwhr? ] 
tan 0 4by*b,3 tan? 0) | A E tr ( 2 tan 6 — by*bp 1/, bwbrh) 


(See Eq. (5) for definition of a.) 


3 
Bex — 3b,2brbw — + 3b,2by*? tan — 2b,*bw tan 0) + ty (=) | 
1 tw [brbw’ barby? 
A +a ( tan 0 — bwbsbr) | + to(*/ + by ) 


by* 
+ by*bp — 
cos 6 


cos 


E (bwhz) + in(*/ bwbr) | 


t(?/sbz? + bw*b, — bwbz?) + tw (/ 3 


Jy 
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tan 6 
M, = — byb;, tan 6) + + — tan 6 — 1/, by’? 
br*by ; 
Jp =t, 5/2 + by? tan? — tan + 2/3 bz? + — | + tw + + ‘ 
cos 
tan? tan 6 by® 
1/3 — 1/2 bpby? ——— + by? be tan? — tan + 1/3 + in, | 
0 cos 6 cos 6 
(2) Stringers with unequal flanges and no lips (see Fig. 5b): 
Jr = ('/s)tw® [bs + bw + dr] 
A = twlba + bw + br] ’ 
t 2 tw? | baby? babwh bywbr? |? 
(See Eq. (5) for definition of a.) 
ction _ lw 4 3 tw? babw? badbwhe || bw? 
R = 6 babw* + 3bsbw*be — | 4 +, + bwhr 
rator 
b 3 t 2 b 2 2 
func- I = tw = + br | 
s) be A + a 2 
@ bei — + — 
(1D) 142 3 4 2 + 4 + + 
M, = (tw/2) baby] 
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Transient Loading of Airfoils at Supersonic 
Speeds 


JOHN W. MILES* 
Unwersity of California at Los Angeles 


SUMMARY 


Using the results for the pressure distribution over a supersonic 
airfoil due to a harmonic velocity distribution and the Fourier 
integrals for step and impulse functions, the pressure distribution 
and lift and moment coefficients have been calculated for the uni- 
form entry of a flat plate into a sharp-edged gust and for the 
sudden deflection of a flap, the latter including the sudden pitch- 
ing of the entire airfoil as a special case. Results are given in the 
form of algebraic equations plus curves for Cr and Cy for M = 1, 
2,4,and ©. 


NOTATION 


Ci = lift coefficient, cf. Eq. (17) 

Cm = mid-chord moment coefficient, cf. Eq. (25) 

= pressure-velocity Green’s function, cf. Eq. (4) 

= cf. Eq. (7) 

= Bessel function of first kind, mth order 

Mach Number 

= free-stream velocity 

= semichord 

= gage pressure 

= 

= time (¢ = 0 when leading edge meets gust or when flap 
is deflected) 

= vertical velocity of airfoil 

= gust velocity, positive up 

= dimensionless chord variable, cf. Fig. 2 

= coordinate of leading edge of flap 

= (M?—1)!4 

= cf. Eq. (5) 

cf. Eq. (8) 

= cf. Eq. (13) 

= cf. Eq. (35) 

= cf. Eq. (9) 

= cf. Eq. (33) 

= x asa variable of integration 


~ 
ll 


air mass density 

cf. Eq. (19) 

cf. Eq. (42) 

= angular frequency of harmonic oscillation 


ll 


VELOCITY DISTRIBUTION 


VELOCITY DISTRIBUTION along an airfoil entering 
a sharp-edged vertical gust at time zero is given by 
the unit step function: 


‘Received October 12, 1947. 
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ing. 
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Fic. 1. Unit step function. 
b 
We 
4 x 


Fic. 2. Velocity distribution on airfoil at (Ut/b) = 1 -—~«. . 


w(x, 4) = wl[Ut — (x + 1)] (1) 


where the airfoil leading and trailing edges are at x = 
+* 1, respectively, U is the relative wind velocity, and 
Wp is the gust velocity. The unit step function is shown 
in Fig. 1. The leading edge of the airfoil meets the 
gust at time zero, and the velocity distribution at Ut = 
1 — xis shown in Fig. 2. 

The unit step function may be expressed analytically 
by the Fourier integral! 


lim 1 
J 


+2 

= 
e+ jw 
(Although Eq. (2) is improper for w = 0, this impro- 
priety may be removed either by expressing Eq. (2) 
as a contour integral or by transforming to Dirichlet’s 
integral,! but neither approach possesses any particular 
advantage in the present analysis.) 
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PRESSURE DISTRIBUTION 


It has been previously shown? that the pressure distribution along a two-dimensional (i.e., infinite span) airfoil 
due to a harmonic velocity distribution w(x, ¢) (airfoil relative to fluid) of angular frequency w is given by 


0G jw [* 

p(x, t) = n+ we, oa NGC, (3) 
Gla, 8) = (4) 

1 = (M/s) (5) 


where the notation of reference 2 has been changed slightly by introducing the variable n. Introducing Eqs. (1) 
and (2) in Eq. (3), partially reversing the order of integration, and taking the limit as e approaches zero yields 


UV? 1 
P(x, t) = [Jo(x, &) — jMT(x, ie} (6) 
In = e J, (nw)dw (7) 
w= 1— [(€ + 1)/UE] (8) 
vy = Mn — ut (9) 


. The integral (7) may be evaluated from the standard result? 
Sore" In = (9? — { — v2)1/2 — jv] /n}” (10) 


where the positive root is taken if n? > v?, and the positive imaginary root if 7? < v?. Using Eq. (10) it is found 
that 


To = 2(q? — v2)-1/2, n? > »? 

= Q, n? < y? (11) 
= —2j(v/n) — > v? 

= 0, n?< 


Substituting Eq. (11) in Eq. (6) yields 


p(x, (x, t) =[1@ (1 (n? — v*)-1/ dé (12) 


+ 1)/Ut] (13) 


where only the real part of the integral is to be included 
in accordance with Eq. (11). Introducing the change 
of variable 


u = v/n = u(x, (14) 


in Eq. (12) and observing that the integral in Eq. (12) 
is real only if « has a magnitude less than unity yields 


P(x, t) = [1(z) — (1 — du] 
(15) 


On carrying out the integration in Eq. (15), it is ex- 
pedient to plot the variation of the variable u as the 
variable is brought from the leading edge (—1) to 
the point at which the pressure is calculated (x). Eq. 
(14) may be regarded as a family of curves for each 
value of x, giving @ as a function of £, where Z@ is the 
family parameter. Actually, it is more convenient to 
plot Ut vs. £, since the curves are then straight lines. 
The points of interest are bounded by £ = —1 and x 
and « = * 1, as shown in Fig. 3. It will be observed 
that the intersections of u = *1 with € = —1 occur 
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+/___Ut=0o 

Ut=0 Ut=2 ! 

Fic. 4. y vs. x as given by Eq. (13). 


at @ = +1/M, respectively while vu = +1 and é = x 
intersect at @ = 0. Hence, for any fixed time, the 
integration is carried out over a horizontal line in the 
shaded area of Fig. 3, and, if @ < 0, the lower limit of 
u in Eq. (15) is —1 and 1(a) = 1, while, if a< 0, the 
lower limit is +1 and 1(z) = 0. The result of integrat- 
ing Eq. (15) is then 


p(x, t) = 0, a< —(1/M) (16a) 

(16b) 

4, z>1/M  (16c) 


Lirt COEFFICIENT DUE To 


Since the chord of the airfoil is 2, the lift coefficient is 
given by (since the pressure jump is —2)): 


oC, Lit) 


Civ = = =~ p(x)dx (17) 


U 2 U 


Introducing the variable @ from Eq. (13), Eq. (17) may 
be written 


(8/4) Cro = (Ut/2)f, plada (18) 
7 = 1 (2/U2) (19) 


The integral in Eq. (18) must be evaluated in three 
parts, corresponding to the three different pressure dis- 
tributions of Eq. (16). To facilitate these integrations, 
Z is plotted vs. x for various Ut in Fig. 4. Substituting 
Eq. (16) in Eq. (18) yields 


4" ™ Mi 


1 
1 
(20a) 
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= (Ut/2) di, ™>(1/M) — (2%) 


The integral involving the are cosine may be evaluated 
by introducing the change of variable 


= (1 — — (21) 


whence 


costudu [ cos“! u 
Mr J-1 (M—w? + 
u(t) 
(1 — (M — aw | 
—— 


(Ut/2) 


(22 
sin (are | ) 


Siu da = [(M — 1)/M] (23) 


where the last integral is given by Peirce 195.4 
Substituting Eqs. (22) and (23) in Eq. (20) yields 


B. _(8\(Ut 

Cw = (2)(2), < M (24a) 
= {00s [ae “(3 


7 > (1/M) 


The values of (8/4)Cz» given by Eq. (24) are plotted 
vs. (Ut/2b) for M’= 1, 2’, 2,4, and © in Fig. 5, where 
bis the semichord. 


(24b) 


= 1, (24c) 


MoMENT COEFFICIENT DUE TO GUST 


The mid-chord moment coefficient (positive stalling) 
is defined by 
M(t) 


oe = = = 
*/2pU?(wo/ U) 2? 


Introducing Eq. (13) yields 
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ft / 


(¥ | 
Fic. 5. (8/4)0C1L/O(wo/U) vs. (Ut/2b) as given by Eq. (24). 


5 Cue = 


Civ — ; (Ut)? (1 — @)p(a)da (26) 


As in the case of Eq. (20), Eq. (26) must be evaluated 
separately for three different ranges of the parameter r. 
The integral involving Eq. (16b) is 


(1 — @) cos“! | di = 
£ cos~' udu —B* O 
mr? J-1 (M — J-1 (M — u)? 
(27) 


where the last integral may then be obtained by dif- 
ferentiating the integral in Eq. (22) with respect to M, 
treating u(r) as independent of MV during this operation. 
Substituting the results in Eq. (26) yields 


Mew = (Zr - 2), 


2 


Ut\?/ 1 u/s 1 


M M 
=0, 17> (1/M) (28c) 


The values of 28Cy. given by Eq. (28) are plotted 
vs. (Ut/2b) for M = 1, 2’, 2, 4, and © in Fig. 6. Al- 
though it might be more logical to plot (8/2) Caw (since 
all supersonic pressure coefficients presumably vary as 
1/8), it is evidently more convenient to use (2/M) as 
the reference quantity. 


TIME TO REACH FULL LiFT 


Full lift and corresponding zero mid-chord moment 
are reached at r = 1/M or 


(U1/28) = =| (29) 
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Fic. 6. 280Cs,/0(wo/U) as given by Eq. (28). 


| M 
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Fic. 7. (Ut/2b) vs. M for attainment of total lift. 
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so that equilibrium is attained as the airfoil just finishes 
entering the gust only for high Mach Numbers, with 
an increasing delay at small Mach Numbers (although 
the linearized theory used in the calculations of course 
breaks down on Mach Numbers too close to unity). 

A plot of Eq. (29) is given in Fig. 7. 


SUDDEN CHANGE IN ANGLE OF ATTACK 


A problem quite analogous to the foregoing is the 
sudden change in angle of attack of some portion of the 
airfoil. If the sudden change of angle occurs over the 
entire chord it corresponds to a sudden gust loading over 


| 
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tion is given by =0,. bh 


The pressure distribution is then given by modifying Eq. (6) to read 


where J) and J, are given by Eq. (11) if only vis replaced by 
vy = Mn (t — b) (33) 
Introducing the change of variable 
w=1-— — §/U(t — b)] (34) 
wo = 1 — + 1)/U(t — (35) 


and assuming that f is independent of & then yields the result 


2 1 
p(x, t) = — ao E — &) — — (1 — w)-! (1 — dw + 


1 
5(w — wo) (1 — M*w?)-! | (36) 


More general results may be obtained by specifying a velocity distribution with an appropriate dependence of h 
on £, the previously treated entry of a sharp-edged gust being a special example. 


SUDDEN FLAP DEFLECTION 


The result, Eq. (36), will be applied to the sudden deflection of a flap extending from x = x to x = 1, in which 
case Eq. (36) becomes 


p(x, 1) = 0, x < (37a) 


p(x, t) = — al 1 - Pm (1 — w)-! (1 — dw + 


RI dw] > x» (37b) 
uo 


= 1 — — x)/Ut] (38) 


The evaluation of Eq. (37) is analogous to the evaluation of Eqs. (15) and (22) and leads to the following results 
(for x > xo): 


P(x, t) = —(pU*/B)a0(B/M), — (1/M) (39a) 


— (pU?/B) a, fo > (1/M) (39c) 
The lift coefficient due to Eq. (36) is found to be 


_ 1 1 
Cras = (1 — (1 — M? diz (1 — M? — 


Bo (1 — as (40) 


the chord, as contrasted to the uniform entry into the w(t, t) = Uaod(é — &) 1 (t — to) (30) 
gust treated in the foregoing. 

The velocity distribution over the airfoil due to a 6(é — &) = 1, E= & 
sudden increase ao of angle of attack ai time fy and loca- (31) 
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= [211 — (B/M), < — (1/M) (41a) 


2(1 — x0) — %0) | (1 — -? |+ 


1 1 
b 
us M 


= 2(1 — x0)/B, < (1/M) (41c) 
= 1 — [(1 — x) /UE] (42) 


It may be observed that, if (8/2)(1 — x0)—'(OC,/0ap) is plotted vs. [Ut/(1 — xo)b] (the time in flap lengths), the re- 
sult is independent of xo. This result is plotted in Fig. 8 for M = 1.2, 2’, and 2. 
The mid-chord pitching moment coefficient due to Eq. (37) is given by 
1 
(1 — M? dj — 


1 
Ho (1 — dito — = ff fio(1 — fo) (1 — (43) 


1/1°— x0?\{/ B 1 (Ut)? 


1 1 — x? 11 (1 |+ 1 (Ut)? _2 


BUt 1 1 
3 1 Ut)(1 — <— (44b 
(1 + Ut)( To”) us” (44b) 
= —(1/8)[(1 — x0?)/2], 7% > (1/M) (44c) 
1O O4 
ra 3 | a. | 
| at lA | AL Neel 
| 26) x:-/ 
|_| 
| 7 2 3 4 5 
= Ue_ | Fic. 9. 280Cy/av as given by Eq. (44) for x = —1. 
| | 
Fic. 8. (8/2)(1 — xo)~! as given by Eq. (41). Ve 
Ox x,=0 
Evidently the behavior of 0Cy,/0 ao is sensitive to the 
value of xo, and it is not possible to draw a universal ‘d 


curve, as in the case of OC,/Oao. In Fig. 9, 28(0Cy + 
Ya») is plotted vs. (Ut/2b) for x». = —1 and for M = a 


1.2, 2", and corresponding to a suddenly imposed 

angle of attack of the entire airfoil. In Fig. 10, », SZ 
— xo?) is plotted vs. (Ut)/(1 — xo)b Ne 

for x = 0 and for M = 1.2, 2", and 2, corresponding 3 — — 
to the sudden deflection of a 50 per cent chord flap. “10 a 


It does not appear that the transient behavior of a super- Fic. 10. 280Cy/Oao as given by Eq. (44) for xo = 0. 
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sonic airplane due to sudden pitching or flap deflection 
will present serious problems unless dynamic sta- 
bility is marginal (except very close to JJ = 1, where 
the foregoing calculations are probably not too reli- 
able). 

The time in flap lengths Ut(1 — x)—'b-! required 
for attainment of full lift and moment is evidently 
given by Eq. (27) and Fig. 7. 


COMPARISON OF GusT ENTRY AND SUDDEN ANGLE OF 
ATTACK 


It is of interest to compare the results of Eqs. (24) and 
(28) with those Eqs. of (41) and (44) for x» = —1, the 
former representing the uniform penetration of a sharp- 
edged gust and the latter representing either the sudden 
application of an angle of attack aj = wo/U or the sud- 
den application of a gust wo over the entire chord at 
time zero. While the eventual lift and moment coeffi- 
cients and time required to attain these coefficients is 
the same in both cases, the manner in which they are 
reached is quite different, the initial lift being finite in 
the latter problem. 


COMPARISON WITH SUBSONIC CASE 


The incompressible transient loading problem dif- 
fers markedly from the supersonic problem because 
of the influence of the wake in the former problem. 
It has been treated by Kiissner,*»* Wagner,’ Bryant 
and Jones,* and von K4rman and Sears.’ Kiissner 
treated the problem both theoretically’ and experi- 
mentally;* Wagner treated the general problem of 
transient loading; Bryant and Jones used Wagner’s 
results to compute the sudden imposition of the gust 
to the entire chord (i.e., the second problem treated in 
the present paper); and von Karman and Sears treated 
the uniform gust entry. 

The results of von Karman and Sears for the lift co- 
efficient are plotted in Fig. 11. Inasmuch as they 
proved that the lift always acts through the quarter 
chord during the gust entry, the quarter chord moment 
coefficient vanishes identically, 


EXTENSION TO MorE GENERAL LOADINGS 


The effects of more general gusts (i.e., not sharp- 
edged and uniform) may be calculated from the fore- 
going results by Duhamel integration.!. A graphical 
method of accomplishing this end has been given by 
Jones.!° 

Bryant and Jones* have included elastic and stability 
considerations in their subsonic treatment, and their 
techniques may be used, together with the present re- 
sults, to include such effects in the supersonic prob- 
lem. 
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Fic. 11. 


(1/2%)0CL/0(wo/U) for subsonic gust entry, 


RECENT WORK 


Since the foregoing was written, an A.A.F. transla- 
tion of a paper by Schwarz!! has appeared which con- 
tains the result, Eq. (16), and a general analysis similar 
to that of reference 2. 
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Plastic Flow of 
Aluminum-Alloy Plate’ 


L. J. KLINGLER? ann G. SACHS} 
Case Institute of Technology 


ABSTRACT 


Commercial hot-rolled 1!/2-in. 24ST aluminum plate possesses 
q marked degree of crystallographic anisotropy. Tensile tests 
at various orientations showed that the yield strength is only 
slightly dependent upon the crystallographic anisotropy. Ratios 
of the transverse strains to the longitudinal strain varied from 
—0.30 te —0.57, while for an isotropic material the ratio would 
be —0.50. 


INTRODUCTION 


WO TYPES OF ANISOTROPY may be encountered in 

annealed metal. ‘“‘Crystallographic anisotropy” 
results from a preferred orientation or “‘texture’’ in the 
crystal structure of the metal, while “mechanical 
anisotropy’ results from a geometrically directional 
arrangement of certain phases, inclusions, porosities, 
cracks, etc. In addition to these two, any cold-worked 
metal develops a third type of anisotropy which may be 
termed ‘‘anelastic anisotropy.’’ This type appears 
to be caused by a certain kind of residual stresses set up 
during the cold working. It can be reduced and seem- 
ingly eliminated by a stress relief anneal. A paper has 
been published recently! which discusses this type of 
anisotropy. 


This paper will deal only with the first of these 
types—namely, the effect of the crystallographic 
ansotropy on the plastic flow properties. 


Crystallographic anisotropy develops gradually in — 


an initially (apparently) isotropic metal in proportion 
to the total reduction. This has been shown repeat- 
edly by means of both tensile tests” * and x-ray investi- 
gations." Annealing a cold-worked metal results 
either in partial retention of the previous anisotropy or 
in the development of a new type of anisotropy. The 
relation between these annealing and cold-work textures 
is not definitely known.‘ 


The degree of anisotropy in annealed metal depends 
upon a variety of factors, the most important of which 
are the following, listed in the order of decreasing im- 
portance: (1) degree of cold work, (2) temperature 
of the final anneal, and (3) temperature of the (penul- 
timate) anneal before the final cold work. 


Received May 13, 1948. 

* The work reported in this paper was made possible by the 
Frederick Cottrell grant received from the Research Corporation, 
New York. 

t Metals Research Laboratory. 


The degree of anisotropy increases as the tempera- 
ture of the final anneal increases and as the temperature 
of the penultimate anneal decreases. Generally, a metal 
reduced by more than 80 per cent in cold rolling or 
drawing, and subsequently annealed, exhibits a pro- 
nounced annealing texture. However, certain metals 
have a low tendency to develop an annealing texture, 
probably because they contain heterogeneous particles 
(impurities, hard phases, oxides, etc.). 

Materials subjected to hot work (rather than cold 
work) may be substantially free from crystallographic 
orientation and anisotropy. However, little is known 
about how hot working affects these phenomena. A 
few investigations on hot-rolled steel strip® 7 show the 
development of preferred orientation to be especially 
evident when the rolling was carried out below the A; 
temperature. The degree of preferred orientation was 
found to increase with the per cent reduction. The 
developed structure was found to be difficult to remove 
by heat-treatment. 

At present, it is not known whether a wrought metal 
product can be obtained with such a low degree of crys- 
tallographic anisotropy that the deviations cannot be 
revealed by the common methods of testing. 

To obtain a uniform piece of metal suitable for ac- 
curate fundamental investigations, it is general prac- 
tice to process it from a cast ingot by a sequence of hot 
working, cold working, and annealing operations. In 
such a wrought metal product the type and magnitude 
of anisotropy, or directionality, is determined primarily 
by the two last operations as discussed above and, to a 
gradually decreasing degree, upon the preceding opera- 
tions. 

The laws of plasticity have been based on the assump- 
tion that a metallic body can be considered as isotropic. 
However, since it is known that metals may be highly 
anisotropic, it is not certain whether the test results 
offered as proof of certain laws of plasticity were in- 
fluenced by anisotropy. Lode’s tests* on tubes under 
internal pressure and longitudinal tension showed de- 
finite evidence of anisotropy, since the radial and cir- 
cumferential strains were not equal for the case of 
pure longitudinal tension (pressure equal to zero). 
Taylor and Quinney,’ in their tests on tubes subjected 
to torsion and tension, discarded materials that devi- 
ated considerably from complete isotropy as measured 
by change in volume of tubes under tension. However, 
all the tubes investigated showed some change in vol- 
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ume, indicating that the radial and circumferential 
strains were not exactly equal. 

No definite information is available on the effect 
of anisotropy on plastic flow. An attempt by Dorn'® |! 
to establish certain relations for sheet possessing dif- 
ferent types of anisotropy has not yet progressed suffi- 
ciently to account for either the anisotropy in yield 
strength or the anisotropy in plastic flow which are 
observed. 


graphic anisotropy of a sheet (or plate) originates in the 
anisotropy of the crystallites. Single crystals are 
highly anisotropic;'? face centered cubic crystals, for 
example, possess flow stresses in uniaxial tension which 
may differ in the ratio of almost two to one, depending 
upon their orientation.'* A single crystal subjected to 
uniaxial tension generally strains so that the strain in 
one transverse direction is nearly zero—i.e., very much 
smaller than in the other transverse direction. There- 
fore, these features of anisotropy should be retained, 
over the limits of normal accuracy, in a commercial 
metal if it is to be considered anisotropic. 

In commercial sheet the directionality in flow stress 
and the unequality of the two transverse strains, 
therefore, can be used as measures of the degree of 
crystallographic anisotropy. 

This paper will report on the effect of the crystallo- 
graphic anisotropy on the plastic flow properties of 
tension specimens taken at various or.entations in a 
commercial 24ST aluminum-alloy plate. 


It must be assumed, at the outset, that crystallo-- 
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MATERIAL AND PROCEDURE 

A distribution of properties was made on a commer. 
cial 11/.-in. 24ST aluminum plate. The working scheq. 
ule for 1'/2-in. plate calls for hot rolling to size from 
the ingot, heat-treating; followed by cold straighten. 
ing. The straightening process, however, amounts to 
less than 1 per cent cold work. 

All tests were made on 1'/2-in. threaded-end tensile 
specimens, as shown in Fig. 1. ‘ 

The orientations of the specimens cut out of the 
plate are shown on the stereographic projection, Fig, 2. 
Points on this projection represent the intersection of 
straight lines from the center of the stereographic 
sphere with the surface of the sphere. The points may 
then be considered to represent the longitudinal axes 
of the specimens. 

The positions of the specimens can be more easily 
identified by passing planes normal to the rolled 
surface of the plate at increments of 15° to the rolling 
direction. Thus, in this system, plane A is defined by 
the rolling direction and the normal direction, plane B 
is defined by a line in the rolling plane at 15° from the 
rolling direction and by the normal direction, and so 
forth to plane G, which is defined by the transverse and 
normal directions. Any specimen can then be located 
by giving the plane designation and the angle of ele- 
vation of the specimen axis from the rolling plane. 

.A number was stamped and a line was scribed on the 
end of each specimen blank after it was sawed from the 
plate, so that a record of its original position in the 
plate could be kept. All specimens were made in dupli- 
cate, the test section in all cases being taken in the 
middle of the plate. Some tests were discarded be- 
cause of faulty technique, and some additional tests 
were made as checks. 

All specimens were reheat-treated after machining, 
in the following manner: solution treatment, 920°F for 
45 min.; cold water quench; aged at least 4 days at 
room temperature. This treatment was carried out to 
standardize the heat-treatment of all the specimens and 
to remove any residual effects of machining. 

The changes in diameter (from which the strains were 
calculated) were measured to the nearest 0.0001 in. 
with a vernier micrometer. 

All of the specimens tested were found to exhibit one 
of the transverse principal strains, either the maximum 
or the minimum change in diameter, in a direction 
parallel to the rolling plane. This strain was desig- 
nated as @, with e, being the strain along the longi- 
tudinal axis of the specimens and ¢; being the third 
principal strain.* 

Stress-strain curves were plotted of stress} versus the 
two principal transverse strains (conventional), as 


* ©, €, and e; are the natural plastic strains, being defined as 
= In(1 +e) ..., ete., where and es are the conventionai 
strains after the elastic strain has been deducted. 

+ All stress values are true stress, defined as load divided by 
actual area. 
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shown in Fig. 3a. The elastic strains were taken out, 
and the stress-strain curves were redrawn with stress 
as a function of the two natural strains, €, and ¢; (Fig. 
3b). A stress-strain curve of stress versus €, was then 
constructed from which the yield strength (0.2 per cent 
plastic strain) was determined, as shown in Fig. 3b. 

A special series of tests was made to check the general 
symmetry of the plate—namely, that the anisotropy 
in each quadrant of the system of coordinates formed 
by the rolling direction, the transverse direction, and 
the normal direction was the same. Only four quad- 
rants needed to be checked, since the test sections were 
all taken from the center plane of the plate. These 
tests were made on specimens located in the four pos- 
sible C planes at an angle of 45° to the rolling plane. 

It can be seen from the stress-strain curves in Fig. 
4 that the belief in the general symmetry: of the plate 
was valid. The values of each property measured were 
equal for all four quadrants, within the limits of ac- 
curacy and scatter. Therefore, only one quadrant, as 
represented in Fig. 2, need be considered. 

Fig. 4 also shows that the degree of scatter en- 
countered, resulting from nonuniformity of the alumi- 
num plate, was small. 


RESULTS AND DISCUSSION 


Yield Strength 


The plate was found to exhibit a definite degree of 
anisotropy, which seems to be typical of aluminum- 
alloy plate. However, all yield strength (0.2 per cent 
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plastic strain) values, Fig. 5, were found to be within 


+5000 Ibs. per sq.in. of an average value. The ob- 
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served directionality in the yield strength was there-. 
fore small, since the scattering due to nonuniformity of 
the plate amounted to approximately +2,000 Ibs. per 
Considering this uncertainty, the yield strength 
values in the rolling plane exhibited a directionality 
similar to that observed in previous investigations for 
severely cold-rolled and subsequently annealed, pure 
aluminum sheet.‘ 
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Strain Ratios 


The strains differed considerably from those required 
for an isotropic metal—namely, that for any orienta- 


tion: 
& = —e,/2 


Curves of stress versus the two 


plane at angles of 0°, 45°, 


transverse strains, 
€, and ¢3, are shown in Fig. 6 for specimens in the rolling 
and 90° to the rolling direc- 
tion and for a specimen in the normal direction. These 
curves show the degree of anisotropy for the principal 
directions in the plate. They also serve to illustrate the 
different types of behavior encountered. The transverse 
strain, ¢s, was the larger transverse strain for most of 
the specimens tested. The longitudinal and transvers 
specimens, Figs. 6a and 6c, followed this pattern. The 
maximum transverse strain for some of the specimens 


was the strain ¢ as illustrated in Fig. 6b for a spect 
men in the rolling plane midway between the long 


tudinal and transverse directions. 


A few specimens 


exhibited almost isotropic behavior in that the two 
transverse strains were practically equal, as shown for 


the normal specimen in Fig. 6d. 
In Fig. 7, 
showing and ¢; as functions of 


the data for the same curves are replotted, 


It can be seen from 


the last graphs that the points generally fall on 4 
the strains encoul- 


straight line. This means that, for 
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tered in the tensile test, the anisotropy remained 
practically constant for the duration of the test—.e., 
independent of ¢1. 

Taking the ratio ¢,/e, as a measure of the anisotropy, 
the observed values are represented in Fig. 8, and again 
in Fig. 9, as one quadrant of a stereographic projec- 
tion. 

The latter representation shows that the values e/e 
do not conform to a high degree of symmetry or to any 
relation that could be expressed by a simple algebraic 
formula. Such a condition is generally encountered 
in the plastic flow of single crystals. Consequently, it 
must be assumed that the texture of the plate possessed 
some preferred orientation. Thus, the structure of a 
plate may be considered as consisting of a mixture of 
crystals, a part of which was randomly oriented with 
the remainder more or less accurately aligned in certain 
preferred orientations. 

The existence of preferred orientation was verified by 
an x-ray diffraction pattern taken in the normal direc- 
tion at the center of the thickness of the plate. The 
material used for the x-ray pattern was reheat-treated 

in the same manner as the tensile specimens. This 
diffraction pattern, Fig. 10, shows preferred orientation 
comparable to that shown by aluminum that has been 
subjected to a considerable degree of cold work.‘ 

It is impossible at present to predict the plastic be- 
havior of a mixture of randomly oriented crystals and 
aligned crystals. Even a comparatively simple struc- 
ture consisting of a symmetrical arrangement of a few 
crystals has not been analyzed in this respect to date. 
Some light is thrown on this problem by previous in- 
vestigations. The largest degree of anisotropy is that 
exhibited by single crystals.'#!5 However, even for 
single crystals, experimentally determined strain ratios 
did not conform to calculated values based upon a 
single slip mechanism, as illustrated in Fig. 11. From 


Fic. 10. X-ray diffraction pattern taken in normal direction 
at center of 1!/2-in. 24ST aluminum plate (integrating diffraction 
camera, Mo. radiation, 2.93 cm. sample to film distance). 
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this data, it is evident that a small amount of imper- 
fection in the crystal lattice will considerably alter the 
strain ratios. 
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Furthermore, the structure of annealed sheet of some 
face centered cubic metals approximates the crystal- 
lographic structure of a single crystal.*!®!7 Such 
sheet is composed of crystals that are aligned with 
their cubic axes nearly parallel to the rolling, trans- 
verse, and normal directions of the sheet. Some authors 
claim that close to 100 per cent cubic alignment (i.e., the 
sheet being almost a perfect crystal) can be obtained 
by suitable processing. 

The anisotropy in plastic flow exhibited by speci- 
mens taken from such sheet in various directions (in 
the plane of the sheet) has been measured by von 
Goeler and Sachs‘ and by Baldwin, Howald, and Ross.'’ 
As expected, these data show that (particularly for 
small strains) the strain ratios vary considerably with 
the orientation of the test specimen as shown in Fig. 
11.* It can be further seen from Fig. 11 that, even 
though the sheet was presumably 100 per cent aligned, 
there is a marked difference in the strain ratios from 
those for single crystals of the same orientation. This 
deviation of the sheet specimens from the theoretical 
behavior of crystals is in the same direction but in a 
larger degree than the deviation shown between theo- 
retical and actual crystals. This may be explained by 
the fact that specimens cut from an aligned sheet 
conform to symmetry more accurately than do speci- 
mens cut from single crystals; hence, more slip sys- 
tems become operative during the deformation. 

In Fig. 8 the strain ratios for copper sheet are com- 
pared with the strain ratios for specimens in various 
directions in the rolling plane of the aluminum plate. 
As to be expected from the X-ray study, this plate evi- 
dently does not display the same type of preferred ori- 
entation as the copper sheet. 

It is obvious, therefore, that both the yield strength 
and the strain ratios are correlated with the crystal- 
lographic anisotropy of the material. 

This investigation’ was conducted on only one lot of 
24ST aluminum plate, so that variations might be ex- 
pected for other lots and for other alloys. 


CONCLUSIONS 


(1) X-ray diffraction showed the presence of a defi- 
nite degree of preferred orientation in aluminum-alloy 
plate. . 

(2) The yield strength values possessed only a slight 
directionality. In the rolling plane the variation of 
yield strength was similar to that exhibited by cold- 
rolled and annealed pure aluminum sheet. 


* However, for large strains all orientations show a tendency 
to approach a strain ratio e./e, = —0.5. This indicates that for 
small strains the sheet acted as a single crystal, operating under 
single slip or double slip mechanism, but when a certain strain 
was reached, different for each orientation, multiple slip began to 
take place. 
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(3) Ratios of the transverse strains to the longi. 
tudinal strains showed a marked degree of directionality 
similar in type to that exhibited by single crystals, 

(4) It is impossible, at present, to predict the plastic 
behavior of a material (variations in yield strength and 
strain ratios) from a known preferred orientation, 
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Full-Scale Free-Flight Ballistic Measurements 
of Guided Missiles 


L. A. DELSASSO,* L. G. BEY,+ anp D. REUYL? 


Ballistic Research Laboratories, Aberdeen Proving Ground 


SUMMARY 


Special problems arising in the instrumentation of full-scale 
free-flight tests are considered, particularly as pertaining to 
guided missiles. Applications of specialized optical and electronic 
equipment are reviewed in detail. 

Under optical instrumentation such problems as the modifica- 
tion and application of acceleration cameras, photo-theodolites, 
and precision-plate cameras are discussed. In addition, the de- 
yelopment of tracking telescopes using refractors and reflectors 
of large effective focal length in conjunction with special mounts 
for studying the behavior of missiles during the entire flight is 
considered; results obtained so far from these instruments are 


“included. Brief reference is made to the methods for de- 


termining ballistic data, by means of high-speed motion-picture 
cameras, for aircraft rockets launched from airplanes. 

Under electronic instrumentation, the use of radar for trajec- 
tory determinations and the modifications necessary on standard 
radar equipment for use in long-range guided missiles are dis- 
cussed. The use of radio-doppler and the application to this field 
are considered in detail, and the problems and advantages en- 
countered in this field are presented. The usefulness of micro- 
wave reflection doppler equipment is also reviewed. Finally, a 
comparison of the results obtained by the various systems of in- 
strumentation is made. 


OPTICAL METHODS 


General 


Y og PROBLEMS to be solved in the measurement of 
guided missiles by optical methods are similar to 
those encountered in the fields of astrometry and plane- 
tary photography. 

The determination of the trajectory of the missile can 
be conveniently made by instruments of moderate size. 
However, for the measurement of missile attitude and 
the observation of special phenomena, large instru- 
ments are needed. This need has led to the develop- 
ment of powerful instruments, similar to astronomical 
telescopes but differing in the design of mounts, drive 
mechanisms, and auxiliary instrumentation. 

In general, the instruments used at White Sands 
Proving Ground do not represent the ideal instrumen- 
tation for the measurement problems to which they 
ate applied. The urgency of the need for data has 
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often been the determining factor in deciding upon the 
modification of existing methods and available instru- 
ments rather than the design of new ones. For obvious 
reasons the principle of overdetermination is applied 
for each set of instruments wherever possible. 


Trajectory Measurements 


For the measurement of position of the missile, as a 
function of time, the following types of instruments of 
moderate focal length were available at the start of the 
A-4 experiments: 

(a) Acceleration or Bowen-Knapp cameras, fixed, 
of high-speed, ribbon frame, motion-picture type. 

(b) Photo-theodolites, tracking, American-made 
Mitchell and Akeley types. 

(c) Cine-theodolites, tracking, German-made As- 
kania type. 

(d) Ballistic cameras, fixed, of single plate, multiple 
exposure type. 

The fixed cameras will, in general, yield results of 
greater accuracy than can be obtained with tracking 
instruments. However, the observable portion of the 
trajectory is limited for a given setup. 

The photo-theodolite is therefore the preferred opti- 
cal instrument for the measurement of trajectories of 
long ranges and high altitudes. 

The Bowen-Knapp camera has been used exten- 
sively for trajectory measurements of rockets and is 
generally useful for problems that require observations 
at a high rate and over a considerable width of field. 
At White Sands two Bowen-Knapp camera systems of 
three stations each are used. For each launching site 
the stations are located approximately 1 mile to the 
South, East, and West, respectively. 

A photograph of the camera is given in Fig. 1, show- 
ing means for adjustment in azimuth and elevation. 
The frame measures 5'/, in. in width, the height being 
variable in stages between 0.9 and 0.15 in. Corre- 
sponding exposure rates are 30 and 180 frames per sec., 
respectively, the lowest rate being used in most cases. 
The frame may be rotated in its own plane by means 
of the cradle mounting to have its long dimension 
roughly parallel to the projected path of the missile. 

The camera may be equipped with a 7- or a 12-in. 
lens and is operated at a fixed exposure time of 0.0001 
sec. 

The timing equipment used with the camera consists 
of a disc rotating at a rate of 10 r.p.s. by means of an 


605 


longi. 
onality 
plastic 
sth and | 
n, 
tion to 
any of 
| 
| 


606 JOURNAL OF THE AERONAUTICAL SCIENCES—OCTOBER, 


Bowen-Knapp camera equipped with timing device, 
mounted in M-14 trailer. 


Fic. 1. 


1,800 r.p.m. synchronous motor operated at a crystal- 
controlled 60-cycle frequency. A portion of the cir- 
cumference of the disc is photographed through the 
camera objective by means of an auxiliary optical 
system. 

The present synchronization of the instruments of a 
set is achieved through the recording of a photoflash 
lamp, operated by the firing signal pulse by means of 
a thyratron circuit. Hence, the accuracy of synchroni- 
zation is determined by the frame speed and amounts 
to sec. 

Since the portion of the trajectory (the first 4,000 
ft.) measured by the Bowen-Knapp system cannot be 
assumed to lie in a plane, the triangulation is carried 
out for individual points by reducing the. measured 
film coordinates to azimuth and elevation angles. 
Measurements are made with a Mann comparator, 
specially built for this purpose, or with a Zeiss tool- 


Fic. 2. Mitchell photo-theodolite. 
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maker's microscope. The probable error of an angular 
measurement is +0.1 mil umtder average conditions, 

The first tracking instruments at White Sands were 
of the Mitchell type (Fig. 2) and are located 6 miles 
South and West of the launching site, respectively, 

The main features of this instrument are a foga] 
length of 12 in., tracking by a single observer, and g 
maximum frame rate of 24 exposures per sec. The dj. 
rection of the optical axis is recorded on each frame by 
means of dial counters operated through gear trains, 
Although possible backlash in the counterdrive gears 
causes this method to be undesirable, a change was 
deemed impractical under the circumstances. 

The instrument has been modified by installing g 
171/2° shutter. By means of a slit in the shutter and q 
lamp and phototube arrangement, a pulse is generated 
at the instant of exposure and recorded on a Brush 
chronograph, together with signals of 1 and 32 pulses 
per sec. The firing signal is recorded on the chrono- 
graph tape also and, in addition, on the film in the same 
manner as for the Bowen-Knapp camera. 

The measurements of the azimuth and elevation cor. 
rections are carried out by a projection method.! The 
space coordinates of the missile are calculated in the 
conventional manner. Calculated velocities are the 
average velocities for the small time intervals between 
measured frames. A detailed description is given ina 
separate report.? The internal consistency of the film 
measurements was found to be of the order of 0.2 mils. 
When taken over time intervals of one sec., velocities 
may be determined with a probable error of about 1 
per cent. The comparison of acceleration measure- 
ments made optically and from telemetered missile 
accelerometer data indicates excellent agreement be- 
tween the two methods.‘ 

The Askania cine-theodolite (Fig. 3) is superior to 
the American instruments because of the method of 
photographing directly the azimuth and elevation 
circles, thereby eliminating errors of the type caused by 
gear backlash in the other instruments. Additional 
features are the larger focal length lenses and a some- 
what larger film frame. When operated in a set, all 
instrument shutters are synchronized by pulses re- 
ceived from a central station. However, a number of 
undesirable features are present—e.g., the inconsistent 
operation of the shutter, the intricate camera mecha- 
nism, and others. 

At the present time one set of Askania instruments is 
operated at stations 7 miles NW, 7 miles NE, and 3 
miles S from the launching position, respectively. 
Two of the instruments are equipped with 60-cm. lenses 
while the third has a Cassegrain reflector of 100-cm. 
focal length. A second set of 60-cm. instruments is 
operated at stations 9 miles NNW, 12 miles NE, and3 
miles S from the launching site. Either a 1, 2, or4 
pulse per sec. signal may be transmitted to the various 
stations for the tripping of the shutters. The shutter 
operation is again recorded in the same manner as was 
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described for the Mitchell photo-theodolite system. A 
timing accuracy of 0.001, sec. is obtained by this 
method. In addition, two instruments are set up for 
special purposes, such as the study of wind velocity and 
direction at high altitudes from vapor or smoke trails 
produced by the missile. Simultaneous recording will, 
in general, be difficult when the angle subtended by the 
trails exceeds the angular field of view. The use of 
wide field cameras is therefore recommended for this 
type of study. 

The Askania instrumentation has been used to obtain 
rapid trajectory data of adequate, though possibly not 
maximum, accuracy. The methods of measurement 
and reduction are similar to those used for the Mitchell 
photo-theodolite records.°® 

In general, the Mitchell instrumentation will cover 
the trajectory of the A-4 through the point of burnout— 
ie., to an altitude of about 30 miles. The Askania 
system will yield data generally to 40 or 50 miles. In 


_ many cases a single instrument may track and record to 


points considerably beyond these limits. The remain- 
ing portion of the trajectory may be determined by cal- 
culation, taking into account the variation of gravity 
and the rotation of the earth.® 

The ballistic camera was developed at the Ballistic 
Research Laboratories for the study of bomb ballistics 
and is capable of giving data of high accuracy. The 
method is useful for other purposes also, provided that 
the light source is of sufficient intensity to register on a 
background that is necessarily dense because of multiple 
exposure. 

In its original form the method called for two instru- 
ments mounted vertically at the end points of a meas- 
ured base line.” The calibration of the instruments is 
carried out by means of stellar positions. 
of the ballistic program are made by means of a system 
of shutters operated synchronously within the set of 
instruments and with the light source—e.g., a high in- 
tensity flashing light in the aircraft. Night programs 
are obviously more desirable, since the observations 
and the calibration may be carried out under almost 
identical conditions, and the shutter system can be 
omitted. 

Further work on the development of this method has 
-included the design of new cameras having larger plates 
(8 by 10 in.) and the modification of methods of meas- 
urement and reduction. 

Initial experiments with a tilted camera method® are 
in progress at White Sands thereby also examining the 
possible use of this method for the study of long-range 
missiles. 

At the present time two stations (Fig. 4) have been 
erected at White Sands at positions 3 miles East and 
West of the launching site. The cameras may be oper- 
ated at rates of 1, 2, or 4 exposures per sec. _Prelimi- 
nary results indicate that the method may be useful by 
photographing the jet and the light emitted by the hot 
Jet rudders. Also, suitable photographs have been 


Observations - 


Fic. 3. Askania cine-theodolite. 


made of vapor trails. It has been difficult up to now 
to obtain satisfactory images of the missile itself. 

A statement of the accuracy obtainable by the tilted 
method would be somewhat premature, in view of the 
early stage af present experiments. However, it may 
be of interest to mention the accuracy of the determina- 
tion of lens distortion, using star positions for the pur- 
pose of calibration. It was found that the dispersion 
was given by an average residual of 2 or 3 sec. of arc. 
Allowing for the effect of other errors of measurement 
it may be possible to achieve an angular accuracy given 
by a probable error of +0.02 mil. Obviously, this 


Fic. 4. Ballistic camera in tilted position, with shutter assembly 
and sky screen. 
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showing recording and tracking instruments. 


value will depend largely upon the accuracy with which 
the image of the missile, or of the jet, may be measured. 


Attitude Measurements 


A study of the attitude of the missile as a function of 
time may be undertaken in a number of different ways. 
The method described here involves the use of ground 
instrumentation exclusively, thereby avoiding problems 
encountered in the telemetering of data measured in 
the missile or the recovery of instruments carried by the 
missile. 

For the study of roll, the only requirement is the se- 
curing of a clear record of suitable markings on the mis- 
sile. However, for the determination of yaw and pitch 
it is evident that as a first approximation the accuracy 
of measurement will be directly proportional to the 
focal length. It is readily seen that the increase of 
focal length will call for a corresponding increase in ob- 
jective diameter, thus causing new problems of engi- 
neering and operation. Furthermore, the image dis- 
turbances caused by the atmosphere tend to become so 
large that ultimately further increase in focal length 
would not yield any higher accuracy of measure- 
ment. 

At the present time there are two instruments of dif- 
ferent experimental type in operation at White Sands. 
The first instrument was designed around available 
optics and mounts. The second instrument was simi- 
larly designed for an existing mount for the purpose of 
testing a different optical system. 

Because of the limitations of this equipment it has 
been difficult to obtain complete records from both sta- 
tions. However, for a number of rounds, aerodynamic 
data during descent were obtained by one instrument 
in addition to measurements of the initial portion of the 
flight. The latter are used to check the functioning of 
the missile controls during burning and to determine 
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Fic. 5. Twin 4.5-in. tracking telescope on modified M-45 mount, 
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the missile oscillations in free, near vacuum flight be. 
yond the point of fuel cut off. 

The major problem encountered in the design of , 
tracking telescope is caused by the combination of re. 
quirements of large focal length and accurate tracking. 

The first telescope that was constructed consisted of 
a 4.5-in. refracting telescope, mounted on an M-45 
50-cal. gun turret.’ In order to obtain a scale came. 
sponding to large focal length, the telescope was equipped 
with an amplifying lens that served to convert the 
60-in. focal length to an effective focal length of 209 
in. This photographic refractor was equipped with a 
35-mm. recording Eyemo motion-picture camera oper- 
ated at a rate of about 15 exposures per sec. A pair of 
20-power binoculars was used for tracking and a visya] 
4-in. telescope of 48-in. focal length was installed to be 
used for visual observations. The latter instrument 
has since been converted into a second photographic 
telescope, also with amplifying lens to give an effective 
focal length of 150 in. 

The instrument was not intended for the measure- 
ment of position, and its counters give approximate 
values only of azimuth and elevation. Reliance js 
placed upon a time record used to tie in with other 
instrumentation for the determination of the position 
of the rocket at the time of each photograph. The 
timing is obtained by means of 60 and | per sec. neon 
lamp pips recorded on the edge of the film. 

_ The instrument (Fig. 5) is located 8 miles SW of the 
launching site. In spite of the small field size it was 
found possible to track missiles of the A-4 type satis- 
factorily. A vital factor in obtaining good images of 
the missile is that of having proper lighting conditions; 


_hence, for a given set of stations the time of firing is of 


the utmost importance. In view of the difficulty of ad- 


hering to a time schedule in experimental launchings, 
it is necessary to operate more than the minimum 
number of two stations to ensure a fair probability of 
good records. 


Fic. 6. Ten-in. telescope on modified-M-45 mount showing mif- 
ror and recording camera. 
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When the construction of the second instrument was 
undertaken, the need for greater light-gathering power 
was satisfied by a parabolic mirror of 10-in. diameter. 
Used as a Newtonian Telescope, this mirror has a focal 
length of SO in., but again an amplifying system was 
added to obtain an effective focal length of 280 in. 
4 35-mm. Akeley motion-picture camera was installed 
as the recording instrument. Photographed on the 
image frame directly are azimuth and elevation dials, 
as well as a time dial. 
motors; the latter, by a synchronous motor operated 
at a crystal-controlled frequency. 

The telescope is mounted on an M-45, 50-cal. gun 
turret, on a 40-mm. gun carriage (Fig. 6), and located 7 
miles NW of the launching site. 

A sample record is given in Fig. 7. Measurements 
of tracking telescope records have been made for a 
number of rounds. In most cases these records have 
been made from individual stations; they therefore 
consist of the measurement of the projection of the 


- missile axis on the plane of the photograph. However, 


adequate synchronization of the two telescopes has 
been accomplished recently, and calculations are in 
progress to determine the directions of the axis of the 
missile in space from the two sets of measurements. 

A preliminary analysis of the results indicates that 
the axis of the A-4 image on the 4.5-in. records may be 
measured with a probable error of =0.7° when the mis- 
sile is at a height of approximately 25 miles. 

Film measurements have been made by means of a 
measuring microscope equipped with a rotating reticle 
or on the projected image.' 


Special Measurements 


The usefulness of a long focus instrument is en- 
hanced by its application to special types of measure- 
ments. These include the ejection of apparatus, the 
separation of rockets in a stage system, the functioning 
of steering rockets, and a number of similar occurrences. 
Furthermore, it would be possible to record light signals 
emitted from the missile as a method of optical tele- 
metering. 

Successful records have been obtained of the tracking 
of apparatus ejected from the missile and brought down 
by parachute. War-head blowoffs, as well as evidence 
of functional failures, have been recorded. 

Photographs showing detailed jet structure are being 
studied in conjunction with the analysis of the func- 
tioning of the missile motor, Spectrographic work on 
the jet has been initiated. 


Future Work 


The need for further research on optical methods of 
measurement and development of instrumentation is 
clearly indicated. 

Improved synchronization will be achieved by the 
installation of pip timing on the Bowen-Knapp cameras 


The former are driven by selsyn ° 


Ten-in. tracking telescope record of A-4 missile made 
from Station N at a slant distance of approximately 8 miles, show- 
ing entire 35-mm. frame with dial readings of azimuth elevation 
and time counters. : 


Fic. 7. 


and the synchronized operation of the Mitchell photo- 
theodolite shutters, both in progress. 

The data recording system of the Askania cine-theo- 
dolites is being modified by the installation of Edgerton 
lamps for illumination of the circles. A project to 
equip a number of these instruments with an optical 
system of large focal length is in the design stage. 

The study of design specifications for a precision 


* photo-theodolite capable of yielding measurements of 


missiles at distances of 100 miles and more is being 
undertaken by the Panel on Test Range Procedures 
and Instrumentation. 

Additional work is carried on to improve shutter 
synchronization of the ballistic camera system. Re- 
garding further research on long focus instrumentation 
and methods of tracking, the construction of a 16-in. 
telescope on a 90-mm. antiaircraft gun mount has been 
completed. This instrument is of Newtonian type 
and will be provided with a means to vary the effective 
focal length between limits of 400 and 1,000 in. Proj- 
ects including the design of a tracking telescope of 
Cassegrain type mounted on a SCR-547 radar set 
pedestal are in the planning stage. A number of addi- 
tional problems have come up in connection with this 
work. They include the development of a beam dis- 
tributing device to permit a single telescope to serve as 
a multiple-purpose recording instrument and the design 
of precision circles to make the instrument more nearly 
self sufficient. Future problems in conjunction with 
still longer ranges will require the development of 
methods of automatic positioning of the instruments 
operated in chain systems. 

With regard to tracking problems, an investigation 
of methods of tracking may lead to the design of special 
optical systems. : 

Of fundamental importance in all projects involving 
photographic recording is the improvement of recording 
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Fic. 8. Block diagram of transceiver. 


techniques. 
mental study of photographic emulsions and light filters 
has been initiated through the courtesy and with the 
cooperation of the Eastman Kodak Company. 

Plans for further research on spectrographic methods 
are being formulated on the basis of successful prelim- 
inary experiments. 


ELECTRONIC METHODS 


Radar Instrumentation . 


In the early part of the war it was discovered that 
many types of radar, and in particular such types as 
SCR-584, were capable of tracking projectiles in flight. 
A great deal of work was done in this field at Aberdeen 
to apply radar to ballistic measurements, When the 
first guided missile testing was undertaken, it was 
natural to extend this work to those tests. As the 
missiles under test increased in range, it soon became 


evident that it was no longer possible to track the | 


missiles throughout their full flight by ordinary reflec- 
tion radar. The solution of this difficulty was thought 
to be the use of beacons that had already been de- 
veloped for a similar use in the case of aircraft. At 
this time the Ordnance Department requested the 
Signal Corps, through the Evans Signal Laboratory, 
to undertake this phase of the instrumentation because 
of the Corps’ wide experience with both the ground 
equipment and the beacons. Since that time the Signal 
Corps has been handling this phase of the instrumenta- 
tion problem at White Sands Proving Ground and has 
been successful in tracking the A-4 missiles throughout 
their entire flight. 

It is expected that the Evans Signal Laboratory will 
report on this phase of the work in the near future and 
only a brief outline will be given here. 

The equipment at the present time installed at White 
Sands consists of two modified SCR-584’s with ex- 
tended range and installed plotting boards. As an aid 
to tracking, especially during the early parts of the 
flight, two M-2 optical trackers have been installed in 
such'a way as to permit them to position the radars on 
the missile during flight. 

The plotting boards give an instantaneous plot of 
the path of the missile with moderate accuracy and 
will be particularly useful to follow the flight of mis- 
siles that are controlled from the ground. 

In addition, the range, azimuth, and elevation angle 
information from each of the radars is recorded photo- 


As a part of a larger project an experi- - 
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graphically on motion picture film which can be ana. 
lyzed later. The information gained from the photo. 
graphic records is sufficient to give satisfactory tra. 
jectory plots but, in general, is not accurate enough to 
give satisfactory first and second differences. 


A DopPLER SYSTEM FOR THE DETERMINATION OF Vp. 
LOCITIES AND POSITIONS OF GUIDED MISssILEs 


General Description 


The beacon doppler system in the most general form 
involves transmitting a radio frequency wave from the 
ground to the missile receiving and amplifying this 
signal, and transmitting back to the ground a second 
signal that bears a constant phase relationship to the 
signal received at themissile. Heterodyning of thetrang- 
mitted and received signals yields the magnitude of the 
radio frequency doppler effect. In this system, instead 
of depending upon energy reflected from the missile 
body, as is the case when using any of the so-called re. 
flection doppler systems, a signal of many times greater 
amplitude is radiated by a beacon in the missile, thus 
greatly extending the range of the system. 

To distinguish between the system to be described 
and other systems that may be developed or used in the 
future, an alphabetical name has been assigned to this 
instrumentation. DOVAP, derived from Doppler 
Velocity and Position instrumentation, will be used 
hereafter in referring to the particular method of apply- 
ing the doppler effect to be described in this paper. 

Basically, DOVAP provides for transmitting a sta- 
bilized continuous radio frequency wave from a fixed 
ground station to the missile, doubling the frequency of 
this transmitted wave in the missile receiver and trans- 
mitting the doubled frequency from the missile to three 
or more fixed ground receivers. Fig. 8 is a block dia- 
gram of a typical missile beacon, which functions asa 
receiver and transmitter in the missile. By amplifying 
and doubling the’ received reference frequency and 
retransmitting this signal, fixed phase relationship is 
obtained and the technical difficulties encountered in 
attempting to receive and transmit on a common fre- 
quency are eliminated. 

The frequency of the missile signal as received at 
each ground station is heterodyned with twice the fre- 
quency of the original transmitted ground signal. 
This is accomplished by receiving this reference fre- 
quency on a separate receiver at each ground station 
and doubling the frequency by the same method used 
in the missile receiver. The beat or doppler frequency 
is produced by mixing the radio frequency output from 
the two receivers. Land wires are used to transmit the 

doppler signals to a common point for the purpose of 
making permanent and continuous records on 35-mm. 
photographic film. Fig. 9 is a block diagram of on 


form of ground station equipment showing the method i 


employed to combine the reference frequency and the 
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retransmitted frequency from the missile to produce 

the doppler beat. 

In Fig. 10 a typical DOVAP field setup is depicted 
with the various ground stations located approximately 
as they are at the White Sands Proving Ground range. 

Each doppler frequency cycle represents a change in 
path length from the transmitter to the missile to the 
ground receiver of one-half wave length of the reference 
frequency. From the differences in the change of path 
length to three or more receivers, the spatial coordi- 
nates of the missile at any instant may be computed by 
solving a group of equations representing the intersec- 
tion of three or more ellipsoids (spheres if the trans- 
mitter is coincident with one of the receivers). 

A detailed analysis of the mathematical concepts of 
the data reduction methods is beyond the scope of this 
report. Such reports have been prepared by T-4 S. 
Kusaka!® and B. Garfinkel.!! Only the generalized 
method will be presented here as it applies to the overall 
problem of obtaining ballistic data for the missile. 

The data obtained by the DOVAP system in their 
final recorded form consist of a long and continuous 
series of doppler cycle sine waves from each receiving 
station and a time calibration of the film record. It 
will be assumed that the sine wave traces are continuous 
from the instant the missile leaves the ground, thus 
establishing one point on the trajectory whose coor- 
dinates are known. Almost any other known point 
would serve as well, and, in practice (since the doppler 
traces are not always continuous from the start of the 
trajectory), a point some few hundred or thousand 
feet along the trajectory, as determined by photo- 
theodolites, is used. 

Starting at the known point, each doppler cycle repre- 
sents a change in the total distance from the transmit- 
ter to the missile and from the missile to one of the 
ground receivers of one wave length at 77 Mc. or ap- 
proximately 12.5 ft. If m is the number of doppler 
cycles corresponding to any point along the trajectory, 
then m\ is the total change of path distance as the 
missile moves from the known point to the point to be 
determined. It is, therefore, necessary to count (from 
the knewn point) the number of doppler cycles in each 
channel corresponding to each desired point along the 
trajectory. It is most convenient to refer to the points 
along the trajectory in terms of the time elapsing from 
the start of the flight. If the value of m is determined 
for successive time intervals—e.g., !/2 sec.—then the 
points determined will give the position of the missile 
at '/,-sec. intervals. 

It may be shown that the coordinates (x, y, 2) of any 
point P on the trajectory may be evaluated by the solu- 
tion of two linear equations and one quadratic equation 
of the form 


x=az+l 
y=be+m 
Az? + 2Bz + Cz = 0 


‘in tabular form. 
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Fic. 9. Block diagram of ground station equipment. 


DOPPLER Sigma, 


F starion. 


Launcnee STATION POT 
I> if 


(OW DOPPLER 
OOPPLER SiGwar 


Fic. 10. 


where A, B, C, a, b, 1, and m are coefficients whose 
values may be determined from known factors involv- 
ing the coordinates of the ground receiving and trans- 
mitting stations, the initial point, and the values of n 
obtained from the film record. The expressions for 
these coefficients are lengthy, and to obtain solutions 
for them by ordinary desk computing machines would 
involve weeks of computer time to arrive at the co- 
ordinates of a 100-point trajectory. 

The solution has been programmed for the ENIAC , 
electronic computing machine, and assurance has been 
given that an 800-point trajectory can be computed in 
the time of flight of the missile or approximately seven 
minutes. Following completion of the problem the 
machine may be arranged to print the coordinate data 
It will further tabulate the first and 
second differences between successive points from which 
the velocity and acceleration may be computed. 


Recording Equipment 

As has been previously pointed out, the DOVAP 
system measures radial change of distance from each of 
the ground stations. An examination of a typical A-4 
trajectory will reveal that, in general, there are three 
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ALTITUDE 


Fimat FUEL CUT OFF. 
FIRST FUEL CUT OFF 


HORIZONTAL RANGE 


DOPPLER FREQUENCY 


or 
RADIAL VELOCITY 
° 


HORIZONTAL RANGE 


Doppler frequency as a function of missile velocity and 
position. 


Fie. 1}. 


points on the trajectory at which the tangent to the 
trajectory is perpendicular to the line from a particular 
ground station to the point on the trajectory. In other 
words, the missile is moving in such a direction that 
the radial change of distance is zero and the doppler 
frequency will also be zero. 

Fig. 11 shows the approximate shape of a typical 
A-4 trajectory. Below the trajectory curve is plotted 
a curve of the doppler frequency throughout the tra- 
jectory. 

The first point at which zero radial change of dis- 
tance occurs is at the launching point or start of the 
trajectory. (For practical purposes it will be assumed 
that all ground receiving stations lie in the same plane 
and that the missile takes off in a direction perpendicu- 
lar to this plane. Actually each ground receiving sta- 
tion may be at a differerit height, above or below the 
launching point, and accurate data reduction depends 
upon precise determination of the coordinates of each 
station.) 

The missile reaches a point just beyond the peak of 
its trajectory before the doppler frequency again be- 
comes zero. 
point marked (C) on the curve the missile is moving 
away from all stations. From (A) to (B) the velocity 
of the missile and, hence, the doppler frequency are 
increasing rapidly because of the thrust of the burning 
jet. At point (B) the fuel is cut off, and the missile 
continues on its upward path with decreasing velocity 
until at last the force of gravity overcomes the vertical 
component of velocity and the peak of the trajectory is 
reached. The missile continues to move away from the 


Between the launching point-(A) and the | 
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receiving stations by virtue of the remaining horizonta| 
component of velocity until the point (C) is reacheg. 
At this point the horizontal velocity and the verticaj 
component of velocity due to the earth’s gravitational 
field combine to give a net radial velocity of zero, 

Between point (C) and the impact point (F) the 
velocity increases steadily, but, because of the geomet. 
ric relation of the trajectory to the receiving stations, 
the radial velocity increases up to the point (D) ang 
then decreases to zero again at the third point, marked 
(E).. Between (E) and (F) the radial velocity rises 
until the missile strikes at the impact point. 

This somewhat detailed explanation of how the radial 
velocity and doppler frequency vary throughout the 
flight of the missile is presented here to show the need 
for d.c. response throughout the recording system. [f 
the response were other than d.c., many cycles of dop- 
pler frequency would fail to record as the frequency ap. 
proaches and becomes zero at each of the points dis. 
cussed above. Furthermore, when the radial velocity 
is low and true velocity is comparatively high, each 
doppler cycle may represent a movement of the missile 
of several hundred or even several thousand feet. This 
is particularly true at the point (£) of Fig..11. Thus 
the loss of a few cycles of doppler, either caused by mis- 
counting of the final record or poor response of the re- 
cording system, may introduce large errors in the com- 
puted trajectory coordinates. After considering other 
parameters that affect the accuracy of the system, a 
permissable error of 1 cycle in 100,000 is allowed under 
present data reduction methods. Accordingly,* the 
recording system has been designed to accomplish 
faithful reproduction of all doppler cycles present at the 
output of the ground receivers. 

There is no known practical method of producing a 
so-called direct inking record of a.c. wave forms at fre- 
quencies as high as 800 or 1,000 cycles per sec. Photo- 
graphic methods are more or less standard for this type 
of recording and were chosen because of their extreme 
flexibility and excellence of results. 

The doppler frequency voltage from each ground 
station receiver is applied to the horizontal deflecting 
plates of a 3-inch cathode-ray tube, four of which are 
mounted side by side on a standard 19-inch rack’ panel. 
Conventional sweep voltages are not used during the 
recording process but may be applied to the vertical 
deflecting plates to spread the image of the doppler 
wave for the purpose of examination and analysis. 
The cathode-ray traces of the four tubes are photo- 
graphed by means of a 35-mim. camera in which the film 
is passed behind the lens at a substantially uniform 
rate. The direction of motion of the image on the film 

is transverse to the direction of motion of the film s0 
that the film motion supplies a time base for the re 
corded wave. In order that this time base may he ac- 
curately calibrated at all points along the record, 4 
voltage of accurately known frequency derived froma 
time standard is caused to trip a small glow discharge 
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tube, mounted on the rack panel between two of the 


cathode-ray tubes, at 10 millisec. intervals. The re- 
curring image on the film of this glow discharge results 
in a continuous row of time dots running parallel to the 
sine wave doppler traces and provides an accurate 
means of timing the events on the record with respect 
to the take-off or zero time. At each half-sec. interval 
following the start of the timing wave, a sharp pulse is 
superimposed on the wave to permit more rapid count- 
ing of long stretches of record. Furthermore, every 
twentieth such pulse is deleted to further facilitate 
counting. These half-sec. pulses are recorded by means 
of another glow discharge recording tube. (Special 
circuits cause the first half-sec. pulse to be transmitted 
from the blockhouse at the instant the A-4 lifts from 
the launching stand. Thisreference mark is considered 
to be zero time for all instrumentations in use at the 
White Sands Proving Ground.) 


The doppler frequency has also been recorded on the 


film in the form of a small dot representing each cycle, 


just as the time standard signal is recorded. Experi- 
ence has proved this type of recording to be less tedious 
from the standpoint of counting the total number of 
cycles in each record, especially when the counting is 
done by hand. With electromechanical methods of 
counting, the dot type of record may also be used suc- 
cessfully, but experience has again led to modifications 
so that the dots are replaced with dashes whose lengths 
are of approximately '/.-cycle duration. The amplified 
doppler signals trigger electronic switch circuits causing 
a glow discharge to take place in Sylvania R1130 re- 
cording tubes, one for each doppler channel. These 
tubes, along with the time standard dot trace tubes, 
are mounted on the cathode-ray tube panel in positions 
to cause their traces to be recorded on the film between 
the various sine wave traces. A representation of a 
typical section’ of a DOVAP recording is given in Fig. 
12. 


Description of Typical Records 


Of the first 20 A-4 launchings at White Sands Proving 
Ground, as of the time of this writing, the DOVAP 
system has been used on the last 17, not having been 
used on rounds 1, 2, or 3. For these 17 rounds the 
DOVAP system has worked comparatively well, and 
reducible records have been obtained on all but three: 
An experimental telemetering transmitter in the missile 
blocked the doppler transmitter on round 5, the doppler 
signals were illegibly weak on round 7, and the missile 
doppler unit was not turned on because of failure of 
control lines to the missile on round 19. 


No record has been perfect, there having been varying 
degrees of noise, interference, and fading in each film. 
Continued efforts of design and maintenance engineers 
have, however, consistently improved the DOVAP 
performance, with attendant lowering of the discontinu- 
Likewise, experience in reading the 


ity of the records. 


ND RK 
Ta SION DOPPLER SIGNAL 
TATION DOPPLER DOTS 


-~STATION DOPPLER SIGNAL 


Fic. 12. 


Typical section of DOVAP film record. 


imperfect film records obtained to date has revealed 
methods of interpolation over record breaks and rough 
spots which make the final results closely approach the 
accuracy that would be obtained from a perfect film 
record. Figs. 13 and 14 show radial velocity and ac- 
celeration curves obtained from the DOVAP film record 
of A-4 round 17, which is the best record obtained to 
date. The mean scatter of the points is approximately 
0.2 cycle per half sec. The probable error of radial 
velocity measurements made with these data is ap- 
proximately 0.3 ft. per sec. at a velocity of 5,000 ft. per 
sec. 

Overall accuracy of the DOVAP system in determin- 
ing spatial coordinates is dependent upon many factors, 
not the least of which is the geometric relationship of 
the observing stations. There are four triplets of sta- 
tions at White Sands Proving Ground (BFK, FKG, 
KGB, and GBF) which may be used to determine the 
trajectory coordinates. From analysis of the geometry 
of the stations it is to be expected that triplets BFK 
and BGK should give good agreement in the deter- 
mination of the x coordinate, while triplets FKG and 
FBG should give good agreement for the y coordinate. 
All triplets should agree in the determination of the z 
coordinate. Below are tabular data for a particular 
A-4 round showing the coordinates of the trajectory 
near its peak as determined by each of the four triplets. 
Examination of these data shows that expectations are 
confirmed by actual measurements. 


Stations 
BKF BGK FKG FBG 
x 65274 65201 66738 64504 
y —17998  —21257 —20224 —20222 
z 612739 612646 612613 612753 


An error in the determination of the distance u, (trans- 
mitter to missile to receiver) for each of the stations 
equivalent to two doppler cycles would account for the 
greatest discrepancy observed above. This distance 
if determined from 


u, = Cy + mdr 


where C; is the distance from the transmitter to the 
launcher to the receiver, m is the number of doppler 
cycles, and ) is the doppler wave length. Thus, survey 
errors in obtaining C, miscounts in m, and improper 
assumptions as to the value of the doppler wave length 
may each contribute to an error in u,. It is estimated 
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that survey errors may reasonably account for the 
equivalent of one-half cycle and errors in assuming the 
yelocity of propagation in the ionosphere may account 
jor about one-tenth of a cycle. Thus, the remaining 
error would appear to be due entirely to a total error 
in the doppler count of 14/1) cycles out of approximately 
100,000 cycles, which seems to be entirely. reasonable, 
considering the many factors that may contribute to 
such a miscount. 
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Letter to 


Sir: 

In the first part of the paper ‘‘On Airfoil Theory and Experi- 
ment” (Alan Pope, July, 1948, JouRNAL OF THE AERONAUTICAL 
Sciences) the author makes no mention of the effect of the 
trailing-edge angle on lift curve slope, and in the last two para- 
graphs of that section it is stated first that: ‘‘The highest lift 
curve slope is noted for the 63000 series, which also has the most 
rearward a.c.,” and, immediately afterwards: ‘In general, it 
can only be said that airfoils whose minimum pressure peaks 
occur well back from the leading edge will have steeper lift curve 
slopes and more rearward aerodynamic center positions than the 
older forward pressure peak airfoils.’’ This seems to be a contra- 
diction, as the 63 series airfoils, although not very old, have for- 
ward suction peaks. 

In the second part of the paper, the section on Drag, the author 
proposes a method of estimating section drag which, it seems to 
me, can be judged only as a rule of thumb on its ability to produce 
the right answer. Of the two major assumptions made, the first 
(that transition occurs at the minimum pressure position) is 
probably satisfactory at high Reynolds Numbers, but, in wind- 
tunnel tests at 2 to 3 million, transition is frequently further back, 
as must be well known to all who have done such tests, especially 
with transition indicators. The method of compensating for the 
initial thickness of the turbulent layer has no theoretical basis and 
can be judged only on results. Pope says that his results ‘“‘came 
within Acg = 0.0008, and usually less than Acg = 0.0004, of the 
value measured in the tunnel,’’ which, for a low-drag section, 
means a maximum error of 20 per cent and a probable error of 
less than 10 per cent, which seems hardly satisfactory. Figs. 1 
and 2 herein show a comparison of the theory with further results 
of tests on 6x; — 012 and 009 sections at R = 34,9 X 108, and 
zero lift, taken from the same report ACR L5CO5 quoted by 
Pope. At 3 million the error is generally worse than 10 per cent, 
but at 9 million (outside the range 1 to 7 million mentioned by 
Pope) the agreement is good. Correction for the fact that the 


006 FIG.1 
4 Theory 
Cy 
004 
x,- O12 > 
x,-009 @ 
002 
03 04 05 
Distonce Of Pressure Minimum From L.E. (x /c) 
R= 3x10 


the Editor 


minimum pressure is at 0.4c, 0.47c, and 0.63c for 64, 65, and 66 
series airfoils, respectively, makes little difference. 

Fig. 3 demonstrates clearly that transition may not coincide 
with the minimum pressure peak, It shows, at R = 3 million 
and 9 million, the transition point on a flat plate having the same 
drag (by Squire and Young’s method) as the experimental drag 
of the sections tested. It will be seen that at 9 million Pope’s 
assumption is justified; at 3 million transition would have to be 
about 0.10c behind the minimum pressure point. 

I suggest that the only conclusion ‘to be reached is that, while 
Pope’s method may give good agreement in certain cases, it 
would be wiser in general to use a method having a better theo- 
retical basis. 


F. CHEERS 
Division of Mechanical Engineering 
National Research Council 
Ottawa, Canada 
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The Optimum Proportions of a Long 
Unstiffened Circular Cylinder in Pure Bending’ 


HERBERT BECKER? 
Fairchild Engine and Airplane Corporation 


SUMMARY 


The problem of determining the radius and wall thickness of a 
minimum weight circular cylinder under pure bending is encoun- 
* tered during the design of a missile fuselage. Selection of any of 
the three common plastic moduli (secant, reduced, and tangent) 
for use in the buckling stress equation will yield approximately 
the same optimum design proportions. Furthermore, small devi- 
ations from the optimum design geometry will lead to minor 
weight increases. A simplified design procedure is presented 
which depends upon the applicability of the secant modulus the- 
ory, for which substantiating evidence is available. However, 
since it is not conclusive, the use of the simplified design pro- 
cedure, which yields a near-optimum design with little effort, 
must be considered uncertain until such a time as further evidence 
is available to support the secant modulus theory. 

In the analysis it is assumed that the buckling stress is linearly 
proportional to ¢/R, in agreement with existing data. Discussion 
of departures from this rule is omitted. 

The development, which is entirely analytic, employs the three- 
parameter method of describing stress-strain curves. This leads 
to expressions for shell radii and weights which are easily com- 
pared for each of the three plastic modulus theories. 


INTRODUCTION 


Ws A LONG cylinder of small radius is subjected 
to a pure bending moment, the wall required to 
resist this load may be quite thick. Therefore, in 
the interest of weight minimization, it is desirable to 
know the optimum proportions of the cylinder cross 
section or, in other words, the combination of shell 
radius and wall thickness which will have a minimum 
area and which will have just sufficient strength to 
withstand the applied moment before the cylinder 
buckles. 

This problem is encountered in the design of missile 
fuselages, which usually are of small radius and are re- 
quired to support large bending moments. The ex- 
pendable nature of missiles, in addition to their per- 
formance requirements, accentuates the weight prob- 
lem considerably. 

Received April 9, 1948. 

* The author wishes to express his appreciation to George 
Parsons, of the Pilotless Plane Division, for the impetus he im- 
parted to the writing of this article by his contention that appre- 
ciable variations in the shape of the material stress-strain curve 
would have little influence upon the optimum design weight. 
Credit is due Ceorge Gerard, of New York University, for origi- 
nating the graphic solution, depicted in Fig. 1, of differential 
equations of the form of Eq. (3) involving stress and modulus as 
the variables. This was done several years ago during discussions 
between Mr. Cerard and the author concerning the problem of 


optimum design. 
¢ Structures Engineer, Pilotless Plane Division. 
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The conditions for an optimum design require an ex. 
pression for the critical buckling stress of a long circular 
cylinder in pure bending. As indicated by the data cop. 
tained in references 1 and 2, a good approximation ty 
this stress is obtained from Eq. (la) for E, = E in the 
elastic range. Consequently, the relationship presented 
in Eq. (la) will be used for the critical buckling stress 
during the development of the optimum design theory, 
An analysis of the applicability of Eq. (1a) is beyond 
the scope of this report. A discussion of possible varj. 
ations in K and the existence of an exponent for t/R 
different from unity are items outside the range of 
this investigation and will be omitted. It will be 
assumed, for the purposes of analysis, that Eq. (la) is 
correct. 

The choice of a plastic modulus cannot be made with 
certainty at this time. However, it is demonstrated 
that, for an optimum design or for small departures from 
the optimum proportions, little weight discrepancy 
exists as a result of varying the modulus between £, 
and EF,, which probably represent the limits within 
which the correct modulus should lie. The choice of 
modulus does have a strong influence upon the applica- 
bility of the simplified design procedure, which was de- 
veloped because of its rapidity and simple data require- 
ments. The material elastic modulus and yield stress 
are the only properties required for the use of the simpli- 
fied design procedure in comparison to a need for the 
complete stress-strain curve (up to, and including, the 
knee) for the more precise theory. 

Evidence for the applicability of the secant modulus 
exists; however, it does not appear to be conclusive. 
Since the simplified design procedure assumes the 
secant modulus to be correct in the plastic range, its 
usefulness is questionable at present. 

The development of the design theory, which is 
analytic as opposed to graphic, utilizes the three-param- 
eter description of stress-strain curves presented it 
reference 3, with the symbols altered to agree with the 
o — € notation for stress and strain. The use of 0/% 
for the material stress-strain curve shape factor occurs 
frequently throughout the paper and is preferred be 
cause it utilizes two physical characteristics of the 
material instead of the abstract » employed throughout 
reference 3. The derivations of the analytic optimum 
stress equations are incidental to the main thesis of the | 
article and, for this reason, have been segregated in the 
Appendix. 
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PROPORTIONS OF CIRCULAR CYLINDER IN 


Ee de, 
when a= 
do 
Oo 7 
Eeo f,(o) 
a 
Fic. 1. Graphic solution of the optimum design equation, 


(dE./do) + (Ee/o) = 0. 


SYMBOLS 


E = elastic modulus, lbs. per sq.in. 

E, = effective modulus (secant, reduced, tangent), Ibs. per 
sq.in. 

K = coefficient in critical buckling stress equation 

Ki = buckling stress coefficient used in reference 2; magni- 
tude isa /E 

Ke = coefficient in optimum radius equation 

L = shell length, in. 


M = external bending moments, in. Ibs. 
R = shell radius, in. ‘ 
W = shell weight, Ibs. 


shape parameter for stress-strain curve 

{=1 + [0.3853/logio (o,/e2)] }[reference 3, Eq. (29)] 
= shell thickness, in. 
= strain, in. per in. 


t 

€ 

p = density, lbs. per cu.in. 

¢ =critical buckling stress = applied bending stress, lbs. per 
sq.in. 

o, = stress, on material stress-strain curve, at = 0.7E 


(Fig. 7), lbs. per sq.in. 

o, = stress, on material stress-strain curve, at E, = 0.85E 
(Fig. 7), Ibs. per sq.in. 

yp = yield stress, as described by the 0.002 set criterion 


Subscripts 
0 = optimum design 
r = reduced modulus theory 
$s = secant modulus theory 


= tangent modulus theory 
Superscripts 


Prime (’) is unprimed quantity multiplied by E,/E. (E, o1, 
and o; are the three parameters used to describe stress-strain 
curves analytically.) 


THEORY 


Fundamental Equation 


In order to determine the optimum proportions of a 
circular shell under bending forces, two conditions need 
to be satisfied: The applied and allowable stresses 
should be equal, and the shell weight must be a mini- 
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mum consistent with the design requirements. These 
conditions may be expressed in the forms 

o = KE,(t/R) (la) 
o = (1b) 
W = 2rpLRt = a minimum (1c) 


The minimization of W results in a differential equa- 
tion in R and ¢ for which there is no solution correspond- 
ing to an optimum design when the stress and modulus 
terms are employed. However, logarithmic differ- 
entiation of Eqs. (1) will yield expressions that can be 
reduced to a differential equation in o and E,. This 
leads indirectly to an optimum design through the prior 
determination of an optimum stress. 

Logarithmic differentiation of Eqs. (1) yields 


(do/o) + (dR/R) — (dE,/E.) — (dt/t) = 0 (2a) 


(do/c) + 2(dR/R) + (dt/t) = 0 (2b) 
dW/W = (dR/R) + (dt/t) = 0 (2c) 
These three equations involve four unknowns. The 


geometric terms may be eliminated, leaving an equation 
involving only o and £,, and their differentials: 


(dE,/do) + (E./c) = 0 (3) 


This is the fundamental equation that must be satis- 
fied in order to obtain an optimum shell design. 

From the stress-strain curve of the shell construction 
material, E, is expressible as E, = fi(o). Since Eq. (3) 
may be written in the form E, = f2(c), the value of o 
which satisfies both these functional relations simul- 
taneously is the optimum operating stress for the 
material. 

Eq. (3) may be solved graphically as in Fig. 1. The 
plot of E, against o is derived from the material stress- 
strain curve with little difficulty. However, it is possi- 
ble to solve Eq. (3) analytically using the three-param- 
eter equations of reference 3. Analytic expressions, 
for the pertinent quantities to be developed below, tend 
to eliminate the errors accrued during cross-plotting. 
In addition, they are useful for evaluating effects caused 
by variations in the quantities appearing in the equa- 
tions. For these reasons the development of the 
optimum design theory will be entirely analytic. The 
symbolism of reference 3 has been altered to conform 
with the o — e notation for stress and strain, and the 
analytic derivations of the optimum stress equations, 
may be found in the Appendix, where equation numbers 
are preceded by capital A. 


Optimum Stress Equations 


The most suitable modulus to use for E, probably 
would be £,, E,, or E,. 

Assumption 1: E, = E,.—Through substitution of 
E, for E, in Eq. (3), a simple optimum-design relation is 
obtainable. Upon application of the stress-strain- 
modulus relationships (E, = o/e, E, = do/de) to the 
following slightly modified form of Eq. (3), 


| 

| 

| 

4, 


dE. Ee dE, de + 1 
do oa 
the solution follows readily. Since 
€ e\de_ € 


(1/e)(E, — E,)(1/E,), + (1/e) = 0 


de de 
from which 


then 


= 2 . (4) 


The graphic determination of o) from Eq. (4) is 
accomplished easily by plotting E,/E, versus o and by 
determining the value of o at which this ratio equals 2. 
The analytic solution of Eq. (3), as presented in the 
Appendix, requires that 


= [7/3(n — (A6) 


The dependence of ¢,,/0; upon the stress-strain curve 
shape parameter, o;/02, is delineated in Fig. 2. For 
reference, Fig. 5 contains a table listing several numer- 
ical values of n for corresponding values of o)/02. 

Assumption 2: E, = E,.—In this case the optimum 
design relationship embracing the various moduli 
[corresponding to Eq. (4) above and Eq. (5) below] 
will be omitted, since it is much too complicated for a 
graphic solution. The analytic solution of Eq. (3) for 
E, = E, leads to the expression, j 


Tr — 9) 
3n(n — 2)? 


(A7) 


This curve is also plotted in Fig. 2. 
Assumption 3: E, = E,.—As is shown in the 
Appendix, 
[((E/E,) — 1] 
By applying this expression to Eq. (3) and including 
the stress-strain-modulus relationships, a correspond- 


ence is obtained among the elastic, secant, and tangent 
moduli, and the shape factor, n, viz.: 


EE,/E,? = 1 — (1/n) (5) 


This equation may be solved graphically in a manner 
analogous to that used for Eq. (4). 


(A9) 


n[(E/E,) — 1] 


The analytic optimum stress equation for E, = E,, 
= [7/3n(n — (A8) 


is plotted in Fig. 2, together with the corresponding 
curves for E, = E,and FE, = E,. 


Comparison Formulas 
Examination of Eqs. (A6), (A7), and (A8) reveals the 
following equalities: 


JOURNAL OF THE AERONAUTICAL SCIENCES—OCTOBER, 


1948 


Fic. 2. Optimum stress as a function of the shape parameters 


(o; and oz) and effective modulus. 
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Fic. 3. Design chart for all three moduli. 


[(2n — 3)/(n — (6b) 


These equations permit rapid evaluation of the 
differences among the magnitudes of oo/a; arising from 
the three choices of £, for a particular material. 


Optimum Design 
Shell Radius.—Eqs. (1) may be combined into a form 
“suitable for the optimum design of a long cylinder in 
pure bending. From Egs. (la) and (1b), 
M/nrR*t = (M/rR*)(R/t) 
and 
R/t = KE,/o 
or 
R? = MKE,/rxo? (7) 
This relation, when modified slightly, provides 4 
simple method for the determination of the shell radius. 
Since 


E,/o? = 


| Sh 


From 


The 

The 
shells 

gentm 
| weight 
| modult 
mum ¥ 


4 
In the . 


From E 


a 


Also, 


i 


much clo 
ratios, 
which de 


618 
| 
| 
| 
= 
or 
| 
3 
( Analyti 
between 
E,= E, 
Accord 
| 


form 
pr in 


es a 
dius. 


then Eq. (7) may be revised to read 
= 


The optimum radius is obtained from Eqs. (7) or (7a) 
when E, = E, and o = oo. Since the right side of Eq. 
(7a) is determined uniquely for any particular value of 
‘ ‘ when o = go, then the optimum design curves of 
} Fig. 3 may be drawn from which it is possible to con- 
struct nomograms for rapid design work. 

The shell thickness is obtainable from Eq. (1b), after 
a has been found through use of the relation 


(7a) 


01(00/01) 


Shell Weight.—According to Eq. (1c), 
W = 2npLRt (1c) 
From Eqs. (1a) and (1b), 

o = KE,(t/R) = M/xR* 


or 
| Rt = (MR/KrE3” 

By substituting R¢ into Eq. (1c) and utilizing Eq. (7), 
4 the shell weight is found to be , 


W = (8) 


The shell weight will be a minimum when o = g». 

The largest weight difference will exist between the 
| shells designed on the bases of secant modulus and tan- 
_ gent modulus theories, with the reduced modulus design 
| weight lying between them. The ratio of the tangent 
| modulus optimum weight to the secant modulus opti- 


mum weight is equal to 
In the Appendix it is shown that 


(9) 


= 1 + (3/7) (A3) 
From Eq. (6a), 
(64/01)"~* = [(n)“ 

or 
= 

Also, 
E 3 ze)" E 


Consequently, E,, = E, and 


(Analytically, E,, ~ E,,. However, since the difference 
between the optimum moduli is less than 1/, per cent for 
vi/o: < 1.200, then for all practical purposes, E,, = 
Ey = En.) 

According to Eq. (8) the optimum-weight ratios 
for the secant and tangent modulus theories will be 
much closer to unity than will the corresponding stress 

ratios. Both these ratios have been plotted in Fig. 4, 
which demonstrates that for oy/o2 < 1.14 (nm > 8) the 
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Fic. 4. Comparison of weight and stress discrepancies between 
the secant and tangent modulus theories. 


maximum weight discrepancy between the two theories 
is less than 10 per cent, whereas the stress discrepancy is 
nearly three times as great. This range of 0/02 
corresponds to that in which most of the high-strength 
metal alloys will be found. 


Nonoptimum Design 


The effect upon W produced by deviations of R from 
the optimum value may be obtained with the aid of 
Eqs. (7) and (8). 


W/Wo = (Euoo/E.o)'" (10) 
and , 
R _ (2: By 
or 
(R/Ro) X (W/Wo) X (¢/o0) = 1 (11) 
The result of plotting W/W, against R/R, is depicted 
in Fig. 5 for several values of 1/2. The curves for the 


secant, reduced, and tangent modulus theories are prac- 
tically coincident for each value of o;/02 and appear as 
single lines in the figure. 

The plotting procedure requires the determination of 
o/o, (for the particular value of o;/o2) and the selection 
of several arbitrary values of o/o,. From these stress 
ratios, E/E) and E/E, (corresponding to oo/o, and 
o/o1, respectively) are computed, using Eqs. (A3), 
(A4), and (A5). Then W/W, may be calculated from 
the expression 


(E/E.)/(o/01) 


| (10a) 
Wo (E/E.4)/(o0/01) 


619 
T 
= ° 
= 
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and R/R, is obtainable from 
R/Ro = } (11a) 


The values of W/W, and R/Rp> calculated in this 
fashion may be plotted as ordinate and abscissa on 
graph paper, as shown in Fig. 5. 

The computations for this plotting procedure yield 
information concerning the dependence of the bending 
stress upon variations in the shell radius. These curves 
have also been plotted upon Fig. 5, from which it is now 
possible to determine the effect of a nonoptimum operat- 
ing stress upon the shell radius and weight. As in the 
case of the weight ratio graphs, all three moduli lead to 
stress ratio curves that lie virtually along one line for 
each value of 


Choice of Modulus 


The effective modulus applicable in the plastic range 
probably is close to one of the three discussed during the 
preceding analysis. Some basis may be found for using 
the secant modulus, in recently published data on the 
buckling stresses developed in axially compressed 
curved magnesium-alloy panels.’ These data show 
good agreement with a buckling stress equation of the 
form 


o = KE,(t/R) (12) 


which is identical to Eq. (1b) with E, = E,. This 
suggests the possibility of applying the secant modulus 
theory to shells in bending, since the shape and size of 
the buckles in an axially compressed cylinder resemble 
closely those formed in a shell of the same geometry 
sustaining bending loads (cf. reference 2). 
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Fic. 5. Weight, stress, and radius variations for all three moduli. 
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Fic. 6. Comparison of test data with secant, reduced, and 
tangent modulus theories for critical buckling stress, 


Additional evidence to support the use of the secant 
modulus may be derived from references 1 and 2, which 
contain data on shells in bending wherein the range of 
buckling strains extends from small elastic values to 
extremely large plastic values. The elastic range 
equation is represented by the straight line of Fig. 6, in 
which the buckling coefficient used in reference 2 js 
plotted as a function of R/t. The symbol K, is used to 
distinguish this coefficient from the K appearing in Eq, 
(la), since K, = K(t/R) in the elastic range. 

The magnitude of K, in the plastic range will have to 
be modified in order to bring the theory into closer 
agreement with the test data. Assuming the secant 
modulus to be applicable, then the modified value of K, 
will be 


= Ki(E,/E) (13) 
The buckling stress equation of reference 2 is 
= (14) 


in the elastic range or upon the straight line of Fig. 6. 
In the plastic range, Kz in Eq. (14) must be replaced by 
K,z’, in which case the new buckling stress equation be- 
comes 


= K,’E = = (14a) 


As is demonstrated in the original paper on the secant 
modulus theory,‘ a’ is the stress taken from the material 
stress-strain curve at the same strain as ¢/E. Thus, 
dividing Eq. (14a) by Eq. (14) leads to a simple expres- 
sion for Kz’: 


= (15) 


Values of K,’ taken from the test data of reference | 
appear in Fig. 6, together with the theoretical curves 
for the secant, reduced, and tangent modulus. 

The stress-strain data for 17S-T aluminum-alloy 
sheet, which appear in reference 6 and which were used 
to construct these curves, do not extend beyond a straill 
40. Consequently, it was 
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necessary to obtain large strain data with the three- 
parameter method [Eq. (Al) in the Appendix]. 

In reference | the test data are presented in the form 
of « plotted against R/t. In reference 2 the values of 
were divided by E (= 107 Ibs. per sq.in.) in order to per- 
mit the plot of A, versus R /t. The same procedure was 
used for the test points of Fig. 6. However, since £ for 
the alloy shown in Fig. 7 is equal to 1.067 X 10’ Ibs. per 
sq.in., the points of reference 1 were divided by this 
value before they were plotted on Fig. 6, which indicates 
the applicability of the secant modulus to be justified 
and bolsters the contention that Eq. (12) is applicable 
to the determination of the critical buckling stress of a 
long cylinder in pure bending. 


Simplified Design Procedures 


A fairly good approximation to the optimum design 
proportions may be obtained by using the expressions 
for Ry which appear on the straight lines of Fig.3. The 
weight discrepancies involved, compared to the mini- 


’ mum value at optimum design, are indicated on the 


graph. The use of the secant modulus theory leads to 
negligible error, as may be easily verified by using Fig. 5. 
On the assumption that o,, = o; (where o,, is the 
stress corresponding to 0.002 set), a further simplifica- 
tion is possible; the yield stress may be used, instead of 
01, in the equations for Ry shown in Fig. 3. In this case 
the expression for the optimum shell radius becomes 


Ro = 
When 


E,.=E, Kr = 0.504 
E,=E, ‘Kr = 0.559 
E,=E, Ker = 0.590 


(16) 


On the hypothesis that the use of the secant modulus 
is permissible, assuming o),) = go will lead to little weight 
error as compared to the minimum design. On the 
other hand, considerable weight discrepancies are realiz- 
able using o,, = oo in conjunction with the reduced 
modulus or with the tangent modulus. This is evident 
from Table 1, which contains a comparison of the weight 
errors for several materials using ¢ = o,, in Eqs. (7) and 
(8) instead of ¢ = o (in which case E, = E,,). Obvi- 
ously, the weight discrepancies are supportable for 
E, = E;. The average secant modulus theory error is 
6 per cent, compared to the minimum reduced modulus 
theory error of 16 per cent. The rest of the errors, ex- 
cept for the 21 per cent /, discrepancy, are 30 per cent 
or more and cannot be read on the graph of Fig. 5. 

A simple, rapid design procedure may be devised on 
the basis of these data. Ry may be computed from 
Eq. (16) using K, = 0.504, or 


Ry = (17) 


The shell thickness then may be determined from 
Eq. (1b) in the form 
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Fic. 7. A stress-strain curve for 17S-T aluminum-alloy sheet. - 


to = yp (18) 


This method has the advantage of eliminating the 
need for a stress-strain curve of the structural material 
to be used for the shell, which is necessary for the de- 
termination of 01, and using the yield stress instead, 
which is one of the fundamental mechanical properties 
usually furnished by the material manufacturer. 


DISCUSSION, 


Simplified Design Procedure 


For design purposes, a rapid, reasonably close 
approximation to the optimum proportions is more de- 
sirable than a precise, slightly more involved calcula- 
tion procedure. Since Eqs. (17) and (18) are readily 
adaptable to the construction of design curves for the 
various structural materials and yield a near-optimum 
shell geometry, it would be worth while to examine, in 
some detail, the factors affecting the applicability of 
these equations. 

It was assumed, in developing the simplified design 
equations [Eqs. (17) and (18)], that o,, = o1 = oo, and 
it was shown that the weight error involved in such an 
assumption is small if = This is apparent from 
Fig. 2, wherein it is demonstrated that o,,/0,; does not 
depart very much from unity and from Eq. (10a), 
which indicates that the error will be reduced consider- 
ably because of the cube root. As is evident from 
Table 1, o,, — o; for the materials listed, and any dis- 
crepancy between the two stresses will be reduced to a 
minute value by the cube root of Eq. (10a). Obvi- 
ously, then, the primary factor controlling the use of the 
simplified design equations is the choice of the secant 
modulus for 

- Although evidence to indicate the accuracy of such a 
choice has been presented, it cannot be considered con- 
clusive because of the nature of the data. The tests 
described in reference 1 were conducted upon both tubes 
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and sheets, for which there undoubtedly was a con- 
siderable variation in the values of 1/02 and oy. 
These errors alone could account for a large portion of 
the scatter in the data. Also, the three theoretical K,’ 
curves were determined from material stress-strain data, 
which, although perhaps typical of 17S-T aluminum- 
alloy sheet, may deviate considerably from the material 
properties of the test specimens. 

Apparently, the solution to the problem of the choice 
of E, can come only through the examination of existing 
test data and perhaps the construction and testing of 
additional specimens. At any rate, it cannot be an- 
swered conclusively now. Consequently, the use of the 
simplified design procedure is not substantiated by the 
test data but only is indicated as a possibility. 


_Limitations of Applicability 


The problem investigated in this paper is the deter- 
mination of the lightest weight shell that will support a 
specific pure bending moment. It was assumed tacitly 
that the pure shell is the lightest type of construction 
which could be found to support the applied moment. 
Within a certain range of moments this will be true. 
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However, as the shell radius is required to increase 
with the applied moment, a critical radius will be 
attained beyond which some other type of construction 
will be lighter. Consequently, it is apparent that only 
one band of the bending moment spectrum has been 
investigated in this report. It is necessary to analyze 
additional types of construction—such as stringer, or 
longeron, stiffened shells—to determine their applicable 
moment ranges. When this has been done, a particular 
bending moment becomes identified with a particular 
type of construction which will yield an optimum design 
for that moment. 

Obviously, the optimum designs obtained for pure 
bending must be modified for shells subjected to addj- 
tional loads from axial, torsional, and transverse shear- 
ing forces. In order to evaluate these effects, it might 
be useful to determine the requirements for an optimum 
design pertinent to each separate loading condition, 
The axial case is important in missile design because of 
the large power-plant thrust developed. However, it 
might be difficult to determine a simplified design pro- 
cedure for axially compressed cylinders, since data in 
the plastic range are difficult to find, as a result of which 


Table 1 = Comparison of weight errors, using 


* Gq for all three plastic modulus theories. 


Aluminum Alloy Sheet Stainless Steel Sheet 

Effective 17s -T8 24,8 = Td 75 $ = Tb Stainless | Type Type 

E (108 psi) 10.67 10.50 9.60 10.42 9.35 28.30 26.00 25.68 
o, (10° pst) 40.10 45.60 35.30 72,60 31.00 191.0 33.10 | 30.6 
G, (10° psi) 34.80 41,25 32.30 69.10 28,25 171.0 30.10 25.20 
yp (10° psi) 41.40 46.00 36.00 71.50 31.75 185.0 37.05 36.70 
1.182 1.106 1.093 1.053 1,088 1.117 1.100 1,218 

n 7.22 9.83 10.95 18.10 1.50 9.08 10.30 5.55 
%/ 9, Reduced +72 +733 +740 787 744 +729 +735 +748 
Tangent 640 663 669 +626 

Secant 35.20 39.75 30.85 65,00 27.10 161.7 25.90 27.85 

(10? psi) Reduced 29,05 33.40 26.10 57.20 23.05 139.2 24.30 | 22,85 
Tangent 25,60 30.70 24.25 54.90 21.50 126.6 22,15 19.15 
Secant 1,175 1.197 1.168 1.100 1.172 1.142 1,282 1,518 
/% Reduced 3.425 1.377 1.380 1,247 1.377 1.330 1,523 1,606 
Tangent 1.45 1.498 1.485 1.279 1.477 1.460 1,672 1,917 

Tangent = = = = 

Secant 1.04 1.04 1.07 1.04 1.07 1.08 1.13 1,07 

W/Wo Reduced 1.19 1.21 - 1.16 - 

a =~ Figure 7 b = Alcoa Aluminum Research Laboratories Data ¢ = Allegheny Ludlum Steel Corporation Data 
d + Smith, R., Wyche, B. H., and Gorr, W. W., = "A Precipitation-Hardening Stainless Steel of the 18 Per Cent Chromium, & Per Cent 
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the applicable modulus would be uncertain. Some 
additional test data would be necessary to clarify this 


point. 
CONCLUSIONS 


The optimum radius ofa long circular shell resisting a 
pure bending moment may be expressed as 


Ry = (16) 
When 
Kr = 0.504 
E,=E, Ker = 0.559 
E.=E, Kr = 0.590 


If the coefficient for the tangent modulus is used in 
Eq. (16) and the secant modulus is found to be appli- 
cable, the weight discrepancy will be 10 per cent or less 
(in excess of the secant modulus theory weight) for most 
structural materials. Evidently, then, knowledge of 
the correct plastic modulus is not necessary for optimum 


_ design. 
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The variation of ¢/o) and W/W, with R/R, is inde- 
pendent of the choice of modulus, and for small de- 
partures from the optimum design geometry the weight 
increase is negligible. Consequently, it is not necessary 
to know which plastic modulus is correct for small de- 
partures from the optimum design. 

The use of the simplified design procedure assumes 
that a) = o,,. This is a reasonable assumption only 
if the secant modulus is applicable in the plastic range. 
The evidence to support this contention is not conclu- 
sive. Consequently, it will be necessary to obtain addi- 
tional support for the applicability of the plastic mod- 
ulus theory before the simplified design procedure can 
be used with complete freedom. Should the applica- 
bility of the secant modulus be demonstrated, then the 
simplified design procedure would provide a rapid, close 
approximation to the optimum shell proportions, re- 
quiring only a knowledge of the applied moment, the 
material elastic modulus, and the material yield stress 
as defined by the 0.002 set criterion. 


Appendix A 


(Derivations of the analytic stress, strain and mod- 
ulus relations are based upon the three-parameter 
method of describing stress-strain curves.) 


Mopvu.Lus RELATIONSHIPS 


From reference 3, Eq. (2), 


= (0/E) + C(c/E)” (Al) 
From reference 3, Eq. (7), 
C(o,/E)"~! = (1 — m)/m, = 3/7 (A2) 


since m, is equal to 0.7 [reference 3, Eq. (22)]. . 
Multiplying Eq. (Al) by E/o leads to 


E/E, = 1+ 


However, 
Consequently, 
E/E, = 1 + (3/7)(¢/a)"~" (A3) 
From Eq. (15) of reference 3, 
E/E, = 1 + (3/7)n(o/01)"~* (A4) 


From Eq. (18) of reference 3, 


(A5) 


OPTIMUM STRESS EQUATIONS 
_ Using the secant modulus in Eq. (3) 
dE, _ —(3/7)(nm — _ —E, 


do [1 + 3/7(¢/0,)"~"]? 
—E/o 
1 + 
or 
(3/7)(n — = 1 + 


from which 
= [7/3(n — 
When E, = E,, 


dE,/do = —4X 
(3/7)n(n — 
[1+ V1 + V1 + 


(A6) 


or 


from which 


‘ 1 -1 
2)?} /(n ) 


= [7(2n — 3)/3n(n — (A7) 


Using E,*= E, in Eq. (3), 
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dE, (3/7)n(m — _ ((E/E,) — 1] = 1) (A9) 
do [1 + 
E/e REFERENCES 
- I + (3/7)n(o/o1)"~! 1 Younger, J. E., Principle of Similitude as Applied to Research 
. on Thin-Sheet Structures, Aeronautical Engineering (Trans, 
or A.S.M.E.), pp. 163-169, October-December, 1933." 


= [7/3n(n — (A8) 


2? Lundquist, E. E., Strength Tests of Thin-Walled Duralumin 
Cylinders in Pure Bending, N.A.C.A. T.N. No. 479, December, 
1933. 

3 Ramberg, W., and Osgood, W. R., Description of Stress-Strain 


INCIDENTAL EQUATION Curves by Three Parameters, N.A.C.A. T.N. No. 902, July, 1943 The 
4 Gerard, G., A Secant Modulus Method for Determining Plate 
The expression to be derived below is used in conjunc- Instability Above the Proportional Limit, Journal of the Aero. sidered 
tion with Assumption 3: E, = £,. nautical Sciences, Vol. 13, No. 1, pp. 38-44, 48, January, 1946, ground 
Revising Eqs. (A3) and (A4) to read 5 Schuette, E. H., Buckling of Curved Sheet in Compression and —. 
Its Relation to the Secant Modulus, Journal of the Aeronautica} cie 
(E/E;) — 1 = (3/7)(¢/01)” ep Sciences, Vol. 15, No. 1, pp. 18-22, January, 1948. simple 
(E/E,) — 1 = (3/7)n(o/01)" 6 Heimerl, G. J., and J. Albert Roy, Column and Plate Com. =. 
ressive Strengths of Aircraft Structural Materials—17S-T Alumi- 
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The Relative Merits of Rotary Compressors 


D. G. SAMARAS* anp R. A. TYLERt 
National Research Council of Canada 


SUMMARY 


The relative merits of various types of rotary compressors are 
assessed in the light of the available test data that could be con- 
sidered sufficiently reliable for the purpose. The theoretical 
grounds on which such an analysis may be legitimately made are 
discussed, and the strict requirements of the theory are relaxed 
sufficiently to euable a graphical comparison to be made on the 
simple basis of small stage efficiency and load factor. The re- 
sulting curves are interpreted with due allowance for the neces- 
sary departure of this representation of the test data from that 
required by the theory. It is concluded that the various types 
of rotary compressors are each characterized by a particular value 
of the load factor at which the efficiepcy for the corresponding 
compressor type is a maximum and departure from which in- 
volves a comparatively rapid loss of efficiency. Furthermore, the 
data analyzed and the theoretical considerations indicate that 
the maximum efficiency increases from one compressor type to 
another of lower characteristic load factor. As a result of this it 
would appear that the axial-flow compressor is, currently, more 
efficient than the centrifugal type to the extent of about 10 per 
cent. 


NOTATION 


- 


= Q/[(x/4)D2(1 — v*)u], flow coefficient 
2gH/u*, pressure coefficient . 

rate of volume flow (cu.ft. per sec.) 

rotor tip diameter (ft.) 

hub-diameter ratio, taken as zero for centrifugal com- 

pressors 

rotor tip speed (ft. per séc.) 

stage head (ft.lbs. per Ib.) 

(2gH)*/*/2+/ aNV/0, load factor 

rotor speed of rotation (r.p.s.) 

Reynolds Number 

Roughness Number 

Mach Number 

circulation (sq.ft. per sec.) 


compressor overall total pressure ratio - 


small stage efficiency 

overall compressor efficiency 
energy added to the fluid (ft.lbs.) 
volume (cu.ft.) 


= Soe 


INTRODUCTION 


i PRESENT CENTURY has seen a shift of emphasis 
from the class of machinery covered by the term 
“reciprocating” to the class that may be referred to as 
“rotary” machinery. Since the turn of the century, 
notable and far-reaching advances have taken place 
in this field, outstanding among which has been the de- 
velopment of the gas-turbine power plant for power and 
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propulsive purposes. At the heart of the gas-turbine 
power plant lies the rotary compressor, fulfilling, as it 
does, a vital pumping function on which the power plant 
depends for its ability to do work. Various types of 
rotary compressors have been developed, the more im- 
portant of which are the axial and centrifugal ‘dynamic 
flow’”’ compressors, which depend for their operation on 
rotating blades or vanes, and rotary displacement com- 
pressors involving a positive displacement of the fluid 
by rotating mating lobes or gear teeth. _ 

The selection of the particular type of compressor 
to be employed in a given power plant is based on an 
engineering compromise between the various demands 
of performance, space, weight, economy, and cost and 
depends on the duty the machine is called upon to per- 
form. Fora particular application one type may be ex- 
pected to meet the requirements of such a compromise 
more satisfactorily than any other. In the field of air- 
craft gas turbines, where the initial demand is for high 
performance at high efficiency, the choice is confined, at 
the present time, to the axial and centrifugal types. 
Each of these compressor types has its own advantages 
and disadvantages. The relative importance of the 
pros and cons is likely, in the future, to undergo con- 
siderable change as the development of new high-speed 
aircraft forms proceeds and advances in turbine design 
and performance continue, but today the choice be- 
tween the two compressor types appears to be.a dif- 
ficult one and a matter for argument. Some engineers 
favor a combination of the two types, and power 
plants embodying this view are in production. The Na- 
tional Research Council of Canada has for some years 
past considered the relative merits of the axial type, and 
its advantages and disadvantages have been discussed 
previously.!}? Much verbal controversy has centered 
around the question of the relative merits of the two 
types, but nowhere in the literature does there appear 
to be a systematic analysis of the problem leading to 
definite conclusions based on a sound theoretical foun- 
dation. The purpose of the present paper is to examine 
the problem in a scientific manner and, as a result, by 
making intelligent statistical use of the reliable experi- 
mental data available, to put the merits of each com- 
pressor in the right perspective. 

For completeness, the examination includes all ma- 
chines that could be categorized as “‘rotary pumps’’ and 
for which reliable data were available. These include _ 
(besides the two types of main interest) turbine pumps, 
Roots and Lysholm blowers, ‘‘angular’’ compressors, 
and propellers. Centrifugal and axial-flow compressors 
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have had, of course, wide application in fields other 
than that of gas-turbine power plants, and these appli- 
cations (e.g., radial superchargers for reciprocating air- 
craft engines, axial-flow cooling fans) have also been 
included. Similar methods to those of the present 
paper have been used previously by Weber’ in examin- 
ing the specific speed of water turbines, centrifugal 
pumps, windmills, and propellers and by Sérensen‘ in 


‘ correlating data on steam turbines. 


SIMILARITY 


The conclusions drawn from a statistical examination 
of experimental data as to the general laws relating the 
behavior of particular machines are legitimate only if 
comparison is made in accordance with the principle of 
similarity. Physical similarity may be said to exist in 
measurable physical phenomena that appear in geo- 
metrically similar constructions and which result from 
the same physical causes. Geometric similarity be- 
tween two constructions implies that each construction 
can be represented by the same construction if appro- 
priate changes are made in the scale of length used 
to measure them. Time similarity can be specified by 
a similar definition. Complete physical similarity re- 
quires, in addition to geometric and time similarity, 
that the dynamic, energetic, chemical, and electrical 
fields of the systems compared also be similar. Com- 
pletely physically similar processes are not only de- 
scribed by the same analytical expressions but are rep- 
resented by the same mathematical functions leading 
to the same characterizing nondimensional parameters. 
In practice, because of the complexity of the processes 
involved, it is necessary to sacrifice a certain amount of 
rigidity in the application of the similarity principle. 
Flow systems are compared which are not strictly simi- 
lar in the complete sense defined above but are similar 
in their essential features. That is to say, only a few of 
the processes involved are considered as important, the 
remainder or secondary processes being disregarded 
with the corrections required to be made by their omis- 


. sion being always kept in mind. Thus, depending 


chiefly on the nature of the problem, the important 
processes are separated from those of secondary im- 
portance and are used as the basis for analysis. For 
the sake of simplicity, it is convenient to deal only with 
the case of a perfect gas; the results obtained will apply 
approximately to saturated vapors. 

The most general characteristics of a rotary compres- 
sor are the flow coefficient yg, the pressure coefficient y, 
and the small stage efficiency 7.. (For the definition 
of these and subsequent symbols see ‘“‘Notation.’’) It 
has been experimentally found, recently, by Habliitzel® 
that, for low Mach Numbers, the performance charac- 
teristics of a rotary compressor are not appreciably af- 
fected by the substitution of water for air as the flow 
medium, and from a nondimensional analysis of the 
flow’ it has been established that 


¥ = fily, Re, M, Pr, Ro) (1) 
No = fly, Re, M, Pr, Ro) (2) 


for systems in which the conditions of essential similar- 
ity, outlined above, obtain. Between various systems 
of such a nature the functions f; and f2 will vary in 
form, although they still will be dependent on the same 
variables. Such systems are considered in this paper 
to be exemplified by the various types of rotary com- 
pressors. 

Considering the flow medium as a particular ideal 
gas, the influence of the Prandtl Number is unimpor- 
tant and may be neglected, although it is important to 
bear in mind that, among the secondary processes ig- 
nored in the above relationships, those involving other 
heat transfer nondimensional parameters (e.g., the 
Nusselt Number Nu And Grashof Number Gr) may as- 
sume importance in flows involving large density and 
temperature variations. Examples involving such 
flows must not be compared with other examples of the 
same system, on the basis of the above relationships, 
unless due allowance be made for this divergence. 

The influence of the Reynolds Number Re is of con- 
siderable importance. A great deal of effort has been 
directed toward deriving the effect of Reynolds Number 
on the characteristics of axial-flow compressors, and a 
large amount of experimental data exists.? Unfor- 
tunately, little work appears to have been done on cen- 
trifugal and other types of compressors in this respect, 
although a fundamentally similar behavior would be 
expected to exist. Fig. 1 shows the results obtained on 
axial-flow compressors. The optimum efficiency in- 
creases rapidly with Reynolds Number up to value for 
Re of about 10°, beyond which the influence of Rey- 
nolds Number, although still apparent, is considerably 
less. The influence of compressibility in these tests was 
avoided by the use of low Mach Numbers, and the 
Reynolds Number used is based on the mean axial 
velocity and the blade chord. 

Mostly because of the difficulty of separating the 
respective influences of the Reynolds and Mach Num- 
bers, quantitative information regarding the latter is 
meager. The results of Encke* show the combined in- 
fluence of the two parameters on compressor characteris- 
tics. Qualitatively, the effect of Mach Number would 
appear to result in a greater sensitivity of the pressure 
coefficient y and the efficiency 7 to a change in flow 
coefficient, with increase of Mach Number. The char- 
acteristic curves of the compressor become steeper and 
the choking mass flow decreases. At high values of the 
Mach Number, greater than a certain critical value de- 
pending on the compressor, the efficiency falls rapidly 
with the large increase of rotor resistance and the dis- 
sipation of energy in shock waves. 

The influence of surface roughness, characterized by 
the nondimensional parameter Ro, has been extensively 
investigated.® !° The roughness number Ro is usually 
taken to be the ratio of the statistical mean height of the 
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surface excrescenses to a principal length—e.g., the 
blade chord. It is well known that a wide range of 
roughness profiles exists in which the various roughness 
profiles have a varying influence on the flow. It was 
found that, at high Reynolds Number, the friction coef- 
ficient for a rough surface is approximately independent 
of the Reynolds Number and depends chiefly on the 
roughness number, being proportionately greater with 
increasing relative roughness. On the other hand, at 
low Reynolds Numbers the friction coefficient is un- 
affected, within limits, by the surface quality and meas- 
ured friction coefficients for surfaces of varying degrees 
of roughness agree closely with the values for a per: 
fectly smooth surface. According to these results, 
therefore, a rough surface may be considered as per- 
fectly smooth, as far as flow resistance is concerned, 
when the roughness elements are contained within the 
laminar part of the boundary layer. Since the thick- 
ness of the laminar boundary layer depends on the 
Reynolds Number, there exists a value of Re for a given 
value of Ro, below which the surface behaves as if it 
were smooth. In this range of Reynolds Numbers, 
therefore, the influence of roughness may be ignored. 
At greater Reynolds Numbers it must be taken into ac- 
count. 

An important standpoint from which the relative 
merits of various compressors may be assessed is that 
of efficiency. Rather than take the overall or the stage 
efficiency as a guide, the small stage efficiency 7. is 
considered in this paper as most suitable for compara- 
tive purposes. The small stage efficiency may be de- 
fined!! as the ratio of the elemental work of compression 
to the elemental shaft input work and corresponds to 
the stage efficiency of a compresssor containing an in- 
finite number of stages and producing the same overall 
pressure ratio as the compressor under consideration. 
It is related to the overall efficiency n1:,, which includes 
irreversible losses through friction, throttling, leakage, 
etc., by the equation 


Not. = [(por/ prs)” D/ We 1] 
(3) 


For centrifugal and axial flow compressors 7. is a few 
per cent higher than the overall efficiency (see Table 1). 

In the past, comparison of diverse groups of rotary 
machinery, such as hydraulic turbines, etc., has been 
made on the basis of specific speed as the measure of 
comparison. In the present examination a factor is 
used, hereinafter termed the Joad factor, which is de- 
fined as the inverse of the specific rotational speed!” 
and given by the expression 


= (4) 


In terms of the nondimensional compressor characteris- 
tics this becomes 


b= — (5) 


For centrifugal compressors the convention is followed 
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of taking v = 0 in the definition of 
= Q/[(r/4)D°(1 — v*)u] 


so that in this case ] = y’’*/y'” 
where 


¢ = Q/[(r/4)D*u] 
Comparing Eggs. (1), (2), and (5) we have 
No = fs(l, Re, M, Ro) (6) 


This equation illustrates the main influences affecting 
the performance characteristics of rotary compressors 
and affords a basis for the statistical analysis of meas- 
ured results. It is, however, only an approximate 
representation of actual conditions, since it neglects to 
take into account the lack of similarity in tip clearance 
losses, boundary-layer formations, etc., which occur in 
practice between compressors of the same type. In 
addition to these effects, a further error is introduced in 
the final analysis by the neglect of similarity in me- 
chanical losses. It is well known that similarity of 
bearings and lubrication is entirely different from the 
hydraulic similarity of compressors and is governed by 
further nondimensional parameters. Because of the 
lack in the available experimental data of the informa- 
tion required to include this effect and because of its 
relatively small influence on the main process of com- 
pression, lack of mechanical similarity has been dis- 
regarded in the results presented in this paper, and a 
small stage efficiency based on an overall efficiency in- 
cluding mechanical losses has been used. _ 

The available data, which are considered sufficiently 
reliable to enable conclusions to be drawn with reason- 
able confidence as to their accuracy, are tabulated in 
Table 1. Where possible, the relevant values of the 
Reynolds and tip Mach Numbers have been included. 
For the centrifugal compressors the Reynolds Number is 
based on impeller tip speed (considered a sufficiently 
close approximation to the actual relative fluid velocity 
at the impeller tip) and the tip diameter. Given an 
infinite number of data for each compressor type, Eq. 
(6) requires, for its adequate representation, an infinite 
series of three-dimensional surfaces for every value of 
the fifth variable. In practice, the data available are 
severely limited in number and have been represented 
two-dimensionally in Fig. 2 by plotting the small stage 
efficiency 7., against the load factor / for each type of 
compressor without discrimination of the variation of 
the other variables in each system. The true variation 
of n.. with / at constant values of the other patameters, 
for each system of common compressor type, is there- 
fore obscured by a certain amount of scatter, the mag- 
nitude of which to a certain extent can be judged from 
the tabled values of these parameters and the foregoing 
remarks concerning their influence. Such judgment or 
reasoned estimation has been exercised in drawing the 
curves individual to each system through the scattered 
points. 
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Until now, the discussion has been concerned only 
with individual systems (namely, particular types of 
rotary compressors) involving a certain degree of simi- 
larity within themselves. It is perhaps relevant here, 
before discussing the implications of Fig. 2, to examine 
briefly the theoretical grounds for differences between 
the separate systems—i.e., to extend the analysis to 
include the comparison of one system (compressor 
type) with another. All types of compressors have 
a common fundamental characteristic—namely, the 
addition of energy to a gas stream. - The efficiency with 
which this process is effected is intimately bound up 
with the compressor type insofar as the flow path and 
attendant opportunities for irreversible loss vary from 
one type to another. 


ENERGY ADDITION TO A FLUID BY A ROTOR 


The first clear presentation of the theory of operation 
of rotary machinery was given by Zeuner'® in 1899. 
This theory is the. well-known hydraulic or filament 
theory postulating a steady and uniform flow through 
the machine. The assumption of steady flow implies 


an independence of time for the fluid flow and the 
energy transfer between rotor and fluid. The condition 
of uniformity implies that every fluid filament enters 
and leaves the machine under similar conditions. 


Thus 


the velocity at entrance is the same for all fluid fila- 
ments— that is, the velocity is uniform in direction and 
magnitude over the plane of entrance. Similar condi- 
tions apply at the exit. These assumptions imply con- 
ditions from which there is, in practice, considerable 
deviation. The theory, however, has been widely 
used in the design of rotary compressors, although re- 
cent advances in aerodynamics and gas dynamics have 
contributed new concepts to its refinement. Keller’ 
investigated the transient and periodic fluctuations 
that invalidate the assumption of steady flow. These 
periodic phenomena have been shown to induce pre- 
mature choking and to have detrimental effects on the 
flow quantity.» The assumption of uniform flow ig- 
nores the velocity and pressure distributions set up by 
the rotor blades and the presence in actual machines of 
extensive blade, hub, and casing boundary layer. © 


Further important contributions to the refinement 
of the original approximate theory have been made by 
Miiller,"*. Spannhake,'* and others, according to 
which the flow in the rotor outside the blades is not 
considered as vortex-free. As in the case of the air- 
plane wing, the blades are considered as discontinui- 
ties in the velocity field and may be replaced by a vortex 
distribution.'” The absolute velocity field may then 
be analyzed into two fields—one uniform and irrota- 
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tional and the other rotational due to the presence of 
the blades. This latter field includes the circulation 
that, according to the wing theory, is responsible for 
the production of forces that add or extract energy to 
or from the fluid. The.energy transfer per unit volume 
between the fluid and the rotor is given by” 


dL /dq = (7) 


Thus, the energy addition or extraction per unit rate of 
mass flow of fluid is equal to the circulation produced 
by the blades. Therefore, the addition or extraction of 
large amounts of energy per unit mass flow of fluid re- 
quires the generation of intense circulation, or, in other 
words, since Eq. (7) can be rewritten 


gH = dI/dt (8) 


the production of large stage heads must be associated 
with high rate of circulation. On the other hand, the 
losses associated with an energy exchange depend on 
the path traversed by the fluid particles and the inten- 
sity of circulation present. A complicated and tortuous 
path of travel involving changes of fluid direction and 
fluid cross section, and the presence of intense circula- 
tion, may be expected to be accompanied by relatively 
high losses. It is reasonable to suppose, therefore, 
that a machine of comparatively high loading will be 
inherently less efficient than a machine producing a 
smaller head per stage, unless the path of flow through 
the latter is considerably more complicated. In prac- 
tice, high stage heads and flow path complication go 
hand in hand. Thus, in the turbine pump, the fluid 
completes a full circle in its travel while suffering a con- 
tinuous series of impacts from successive vanes with the 
consequent generation of intense circulation. The de- 
livery head is high (see Table 1). In the centrifugal 
compressor the loading is lower, and the fluid path, al- 
though involving at least two 90° bends and consider- 
able change in flow area, is less complicated. The 
fluid particles are acted upon by the blades only in the 
relatively short distance of their travel from the im- 
peller eye to the impeller tip, with a consequently 
smaller generation of circulation and vorticity. In the 
axial-flow compressor with a comparatively small 
energy transfer per stage the circulation is lower still, 
and, furthermore, the flow path is comparatively direct 
and uncomplicated. 

There are grounds for expecting therefore that, in 
practice the turbine pump will be inherently and in- 
evitably less efficient than the centrifugal compressor, 
which in its turn will be similarly less efficient than the 
axial-flow compressor. 

Elaborating this statement, it is reasonable to sup- 
pose that there exists an upper limit to the efficiency of 
a particular type of compressor which is intimately 
bound up with the nature of that compressor. 

This efficiency represents an ideal efficiency, greater 
values than which are impossible to achieve in practice 
unless the particular compressor type with which it is 
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associated is abandoned. Furthermore, it is to be ex. 
pected that this ideal efficiency will decrease from one 
compressor type to another of greater stage loading, 

The final appeal, however, must be made from theory 
to experimental fact, and we must turn to recorded, 
measured results, employing the theoretical indications 
only to help us interpret these results wisely. 


DISCUSSION OF RESULTS 


As previously mentioned Fig. 2 shows the small stage 
efficiency 7.. plotted against the load factor / for the 
available experimental data covering various types of 
compressors. Curves have been drawn through the 


scattered points corresponding to each compressor type. 


on a basis already discussed, involving recognition of 
the varying values of Re and M in each group. It isat 
once apparent that the experimental points tend to be- 
have in accordance with the simple, theoretical con- 
siderations of this paper to a remarkably close degree of 
approximation. The points pertaining to particular 
types of compressors fall into distinctly separate 
groups, each characterized by a definite, reasonably 
narrow, band of load factor. In each group, as re- 
quired by Eq. (6), the variation of efficiency with load 
factor is similar, and, for each group, a particular value 
of load factor exists for which the efficiency is maximum. 
These ‘‘optimum’”’ load factors range from about 40 
for turbine pumps, through about 20 for rotary displace- 
ment blowers, 5 for centrifugal compressors, to about 1 
for axial flow compressors. Also, in accordance with 
the indications of the previous section, the maximum 
efficiency increases from one compressor type to an- 
other of lower characteristic load factor. In view of 
this excellent correlation with theoretical considerations, 
it is considered that the practical ideal efficiency corre- 
sponding to a given compressor type, the existence of 
which was mooted in the previous section, may in fact 
be a real concept, and the authors have ventured to 
include in Fig. 2 an envelope curve representing their 
considered opinion, based on the experimental evi- 
dence, of the variation of this ideal efficiency between 
various types of compressor. It is at once obvious, if 
this curve is accepted, that the axial-type compressor 
offers the greatest possible returns for development and 
appears to have a margin of efficiency in hand over its 
nearest rival, the centrifugal compressor, of approxi- 
mately 10 per cent. 

For completeness, points have been included to rep- 
resent the propeller. For this group the efficiency vati- 
ation with load factor is obscure. This efficiency, for 
the sake of a consistent comparison, has been taken to 
be the aerodynamic efficiency of the propeller under 
static conditions, given by 


Naero. = V 2/4 (Cr’*/ Cp) (9) 


where Cr and Cp are the usual thrust and power coef- 
ficients, measured under static conditions. Although 
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RELATIVE MERITS OF ROTARY COMPRESSORS 


the envelope curve has been extended to include these 
pints, it is perhaps doubtful whether they sepresent a 
gecific type of rotary compressor distinct from the 
wial-flow type. It is interesting to note, however, 
that they are characterized by a definite load factor, 
ghich is of the order of 0.3. 

The position of the point representing the Heinkel- 
Hirth-Oll angular compressor is interesting in that, as 
yould be expected, it falls between the centrifugal- and 
axial-type groups with a load factor of 2.8. Asa point 
of interest, the load factor of the D.V.L. supersonic 
cmpressor'® was computed to be 2.3. It was not in- 
duded in Fig. 2, however, as being unrepresentative 
because of the very low efficiency resulting from funda- 
mental errors in its design. 


CONCLUSIONS 


Every type of rotary compressor is characterized 
by a distinct, relatively narrow, band of load factor. 

For each type of rotary compressor a particular load 
factor exists for which the efficiency is a maximum. 

The maximum efficiency is intimately linked with the 
stature of the particular type of compressor and in- 
ceases from one type of compressor to another of 
lower characteristic load factor. 

On the basis of recorded, measured results, the axial- 
fow compressor would appear to be, currently, more 
dicient than the centrifugal type, to the extent of 
about 10 per cent. 
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ized reference books to the preparation of exhaustive 
bibliographies, digesting of reports, and general surveys of 
any aeronautical subject. Some of the available services 
are: 


| 

Bibliographies on any aeronautical subject. 
Reports on any aeronautical subject. 

Digests of aeronautical books, papers, periodicals, and refer- | 

ences. | 

Translations. | 

| 

| 


Engineering investigations of special aeronautical subjects. 

Biographies of individuals engaged in aeronautics. 

Photostats of any acronautical or general enginecring 
material. ‘ 

Microfilms made on special order. 

Photographs made from the Institute's photographic collec- } 


tion. 
Drawings and tracings made. 


which comes within the scope of the Service will be ac 
cepted. Special arrangements may be made for work re- 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any material in the collections of 
the Institute may be ordered at standard photostat rates. 


In addition to the services mentioned any commission 
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